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Abstract. Integrable theory is formulated for correlation functions of character- 
istic polynomials associated with invariant non-Gaussian ensembles of Hermitean 
random matrices. By embedding the correlation functions of interest into a more 
general theory of t functions, we (i) identify a zoo of hierarchical relations satisfied 
by T functions in an abstract infinite-dimensional space, and (ii) present a technol- 
ogy to translate these relations into hierarchically structured nonlinear differential 
equations describing the correlation functions of characteristic polynomials in the 
physical, spectral space. Implications of this formalism for fermionic, bosonic, and 
supersymmetric variations of zero-dimensional replica field theories are discussed 
at length. A particular emphasis is placed on the phenomenon of fermionic-bosonic 
factorisation of random-matrix-theory correlation functions. 
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1. Introduction 

1.1. Motivation and definitions 

Correlation functions of characteristic polynomials (CFCP) appear in various fields of 
mathematical and theoretical physics, (i) In quantum chaology, CFCP (i.a) provide 
a convenient way to describe the universal features of spectral statistics of a particle 
confined in a finite system exhibiting chaotic classical dynamics (Bohigas, Giannoni 
and Schmit 1984; Andreev, Agam, Simons and Altshuler 1996; Miiller, Heusler, Braun, 
Haake and Altland 2004) and (i.b) facilitate calculations of a variety of important 
distribution functions whose generating functions may often be expressed in terms of 
CFCP (see, e.g., Andreev and Simons 1995). (ii) In the random matrix theory approach 
to quantum chromodynamics, CFCP allow to probe various QCD partition functions 
(see, e.g., Verbaarschot 2010). (iii) In the number theory, CFCP have been successfully 
used to model behaviour of the Riemann zeta function along the critical line (Keating 
and Snaith 2000a, 2000b; Hughes, Keating and O'Connell 2000). (iv) Recently, CFCP 
surfaced in the studies of random energy landscapes (Fyodorov 2004). (v) For the 
role played by CFCP in the algebraic geometry, the reader is referred to the paper by 
Brezin and Hikami (2008) and references therein. 

In what follows, we adopt a formal setup which turns an n x n Hermitian matrix 
?C = into a central object of our study. For a fixed matrix ?C, the characteristic 
polynomial det„ (cr — Jt) contains complete information about the matrix spectrum. To 
study the statistics of spectral fluctuations in an ensemble of random matrices, it is 
convenient to introduce the correlation function n„|p(<j; k) of characteristic polynomials 

n„|p(^;/*) = /ndet«°(ca -?<:)) . (1.1) 

Here, the vectors <; = (ci, • • • , Cp) and k = (ki, • • • , Kp) accommodate the energy and 
the "replica" parameters, respectively. The angular brackets (/(Jt))^^ stand for the 
ensemble average 

{fm}:K = j df^nmm) (1.2) 

with respect to a proper probability measure 

dUni^) = Pn{n) (D„J€), (1.3) 

n n 

{'Dn'K) = Y[ d'Kjj Y[ c^R-eJfjfe (flmJC^-fc (1.4) 

j=i j<k 

normalised to unity. Throughout the paper, the probability density function P„(?£) is 
assumed to follow the trace-like law 

P„(JC) = exp [-tr„ y(JC)] (1.5) 

with y(5C) to be referred to as the confinement potential. 

There exist two canonical ways to relate the spectral statistics of JC encoded into 
the average p-point Green hmction 

<^n|pW = (ritr„(fe-?C)-M (1.6) 

to the correlation function H„|p(^; k) of characteristic polynomials. 
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• The supersymmetry-like prescription (Efetov 1983, Verbaarschot, Weidenmiiller 
and Zirnbauer 1985, Guhr 1991), 

Gni,M = (n ^) <t(^' ^'), (1.7) 

makes use of the correlation function 

^ni,t' = ( n detn(Ca - ?C) fl ^et-'i^ - ?€)) (1.8) 

\a=l /3=1 / ^ 

obtainable from Tln\q+q' (^, k, k') by setting the rephca parameters k and k' to 
the integers ±1. 

• On the contrary, the rcplica-Uke prescription (Hardy, Littlewood and Polya 1934, 
Edwards and Anderson 1975), 

Gn\pi<i) = (^U lim^K-'-^^ n„|p(^;K), (1.9) 

entirely relies on the behaviour of the correlation function n„|p(<^; k) for real-valued 
replica parameters, k G MP, as suggested by the limiting procedure in Eq. (1.9). 
In this case, the notation n„|p(^; n) should be imdcrstood as the principal value 
of the r.h.s. in Eq. (1.1). Existence of the CFCP is guaranteed by a proper choice 
of imaginary parts of <;. 



A nonperturbativc calculation of the correlation function n„|p(^;/t) of characteristic 
polynomials is a nontrivial problem. So far, the solutions reported by several 
groups have always reduced n„|p(^;K) to a determinant form. Its simplest - Hankel 
determinant - version follows from the eigenvalue representation | of Eq. (1.1) by virtue 
of the Andreief-de Bruijn formula [Eq. (3.3) below] 



n„|p(<j; k) = n! ^ det„ 



Jk 1 



(1.10) 

a=l J 0<j,fc<n-l 

Here, V„ denotes a volume of the unitary group U{n) as defined by Eq. (1.17). 
Unfortunately, the Hankel determinant Eq. (1.10) is difficult to handle in the physically 
interesting thermodynamic limit: finding its asymptotics in the domain n 3> 1 remains 
to a large extent an open problem (Basor, Chen and Widom 2001, Garoni 2005, 
Krasovsky 2007, Its and Krasovsky 2008) especially as the integral in Eq. (1.10) has 
unbounded support. 

For K integers, k ^ IP , so-called duality relations (see, e.g., Brczin and Hikami 
2000, Mehta and Normand 2001, Desrosiers 2009 and references therein) make it 
possible to identify a more convenient determinant representation of n„|p(^; k): Apart 
from being expressed through a determinant of a reduced size (see below), such an 
alternative representation of CFCP displays an explicit n-dependence hereby making 
an asymptotic large-n analysis more viable. For instance, the correlation function 

nH,+,' is, m, m') = ( n deC" (?„ - ?€) fl deU'"'' (?^ - ?£)\ (1.11) 

\a=l /3=1 I ^ 



\ See Eq. (2.1) and a brief discussion around it. 



1 Introduction 



5 



with m e and m' G can be deduced from the result f by Fyodorov and Strahov 
(2003) 



[^"-9'+'=(S')]j = l,...,g; fc=0,...,,+g'-l 
[''■Tl-g'+fe(s)]j = l,...,g'; fc=0, ■■■,«+<?'-! 



(1.12) 



by inducing a proper degeneracy of energy variables. The validity of this alternative 
representation (which still possesses a determinant form) is restricted to q' < (Baik, 
Delft and Strahov 2003). Here, 

A,(^) = det, [.t'] = n - ^a) (1-13) 

a</3 

is the Vandermonde determinant; the two sets of functions, 7rfe(<r) and hk{<;), are the 
average characteristic polynomial 

7rfe(0 = (detfe(^--?£))j, (1.14) 
and, up to a prefactor, the average inverse characteristic polynomial | 

hk-i{<i) = Cfe.i (detfe !(<? - IK))^ . (1.15) 
Finally, the constant Cn,g' is 

ik'/21-L«'/2J (n-g')!V„_,, N„ ' ^ ^ ^ 

where V„ is a volume of the unitary group U(n), 

n(n-l)/2 

V. = ^^. (m, 

The result Eq. (1.12) is quite surprising since it expresses the higher-order spectral 
correlation functions G„|p(<?) in terms of one-point averages (Gronqvist, Guhr and 
Kohlcr 2004). 

For K reals, k G W, the duality relations are sadly unavailable; consequently, 
determinant representation Eq. (1.12) and determinant representations of the same 
ilk (see, e.g., Strahov and Fyodorov 2003, Baik, Deift and Strahov 2003, Borodin and 
Strahov 2005, Borodin, Olshanski and Strahov 2006, and Guhr 2006) no longer exist. 

The natural question to ask is what structures come instead of determinants? This 
question is the core issue of the present paper in which we develop a completely 
different way of treating of CFCP. Heavily influenced by a series of remarkable works 
by Adler, van Moerbeke and collaborators (Adler, Shiota and van Moerbeke 1995, 
Adler and van Moerbeke 2001, and reference therein), we make use of the ideas of 
integrability § to develop an integrable theory of CFCP whose main outcome is an 
implicit characterisation of CFCP in terms of solutions to certain nonlinear differential 
equations. 

As will be argued later, such a theory is of utmost importance for advancing 
the idea of exact integrability of zero-dimensional replica field theories (Kanzieper 

t See also much earlier works by Uvarov (1959, 1969). Alternative representations for n'^^^J^^, {<;,<;') 
have been obtained by Strahov and Fyodorov (2003), Baik, Deift and Strahov (2003), Borodin and 
Strahov (2005), Borodin, Olshanski and Strahov (2006), and Guhr (2006). 

I Making use of the Heine formula (Heine 1878, Szego 1939), it can be shown that 7rj,(c;) is a monic 
polynomial orthogonal on R with respect to the measure dp,{<;) = exp[— y(i;)] d?. The function hk{<;) 
is its Cauchy-Hilbert transform (see, e.g., Fyodorov and Straiiov 2003): 



rfA('?') 

For a review on integrability and matrix models, the reader is referred to Morozov (1994). 
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2002, Splittorff and Verbaarschot 2003, Kanzieper 2005, Osipov and Kanzieper 2007, 
Kanzicpcr 2010). In fact, it is this particular application of our integrable theory of 
CFCP that motivated the present study. 

1.2. Main results at a glance 

This work, consisting of two parts, puts both the ideology and technology in the first 
place. Consequently, its main outcome is not a single explicit formula (or a set of them) 
for the correlation function n„jj,(<i; k) of characteristic polynomials but 

• a regular formalism tailor-made for a nonperturbative description of n„|p(^;K;) 
considered at real valued replica parameters k gMP, and 

• a comparative analysis of three alternative versions of the replica method 
(fermionic, bosonic, and supersymmetric) which sheds new light on the 
phenomenon of fermionic-bosonic factorisation || of quantum correlation functions. 

More specifically, in the first part of the paper (comprised of Sections 2, 3 and 4 written 
in a tutorial manner) we show that the correlation function n„|p(^; k) of characteristic 
polynomials satisfies an infinite set of nonlinear differential hierarchically structured 
relations. Although these hierarchical relations do not supply explicit (determinant) 
expressions for n„|p(^;K) as predicted by classical theories (which routinely assume 
K G Z^), they do provide an implicit nonperturbative characterisation of n„|p(^;/«) 
which turns out to be much beneficial for an in-depth analysis of the mathematical 
foundations of zero-dimensional replica field theories arising in the random-matrix- 
theory context (Verbaarschot and Zirnbauer 1985). 

Such an analysis is performed in the second part of the paper (Section 5) which 
turns the fermionic-bosonic factorisation of spectral correlation functions into its 
central motif. In brief, focussing on the finite- average density of eigenlevels in the 
paradigmatic Gaussian Unitary Ensemble (GUE), we have used the integrable theory 
of CFCP (developed in the first part of the paper) in conjunction with the Hamiltonian 
theory of Painleve transcendents (Noumi 2004) to associate fictitious Hamiltonian 
systems Hi {P{t),Q{t),t} and {P{f),Q{t),t} with fermionic and bosonic replica 
field theories, respectively. Using this language, we demonstrate that a proper 
replica limit yields the average density of eigenlevels in an anticipated factorised form. 
Depending on the nature (fermionic or bosonic) of the replica limit, the compact and 
noncompact contributions can be assigned to a derivative of the canonical "coordinate" 
and canonical "momentum" of the corresponding Hamiltonian system. Hence, the 
appearance of a noncompact (bosonic) contribution in the fermionic replica limit is no 
longer a "mystery" (Splittorff and Verbaarschot 2003). 

1.3. Outline of the paper 

To help a physics oriented reader navigate through an involved integrable theory of 
CFCP, in Section 2 we outline a general structure of the theory. Along with introducing 
the notation to be used throughout the paper, we list three major ingredients of the 
theory - the r function Tn^ (<?, k; t) assigned to the correlation function n„|p(^; k), the 
bilinear identity in an integral form, and the Virasoro constraints and further discuss 
an interrelation between them and the original correlation function n„|p(^;K.). Two 
integrable hierarchies playing a central role in our theory - the Kadomtsev-Petviashvili 
and the Toda Lattice hierarchies for the r function - are presented in the so-called 

II A quantum correlation function is said to possess the factorisation property if it can be expressed in 
terms of a single fermionic and a single bosonic partition function (Splittorff and Verbaarschot 2003, 
Splittorff and Verbaarschot 2004). 
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Hirota form. Both the Hirota derivative ^ and Schur functions appearing in the above 
integrable hierarchies arc defined. 

Having explained a general structure of the theory, we start its detailed exposition in 
Section 3. In Section 3.1, a determinant structure of the r function is established and 
associated matrix of moments and a symmetry dictated scalar product are introduced. 
The bilinear identity in an integral form which governs the behavior of r function is 
derived in Section 3.2. The bilinear identity in Hirota form is derived + in Section 
3.3. In Section 3.4, the bilinear identity is "deciphered" to produce a zoo of bilinear 
integrable hierarchies satisfied by the r function; their complete classification is given by 
Eqs. (3.41) - (3.47). The two most important integrable hierarchies - the Kadomtsev- 
Petviashvili (KP) and the Toda Lattice (TL) hierarchies are discussed in Section 3.5, 
where explicit formulae are given for the two first nontrivial members of KP and TL 
hierachies. Section 3.6 contains a detailed derivation of the Virasoro constraints, the 
last ingredient of the integrable theory of characteristic polynomials. 

Section 4 shows how the properties of the r function studied in Section 3 can be 
translated to those of the correlation function H„|p(<j; k) of characteristic polynomials. 
This is done for Gaussian Unitary Ensemble (GUE) and Laguerre Unitary Ensemble 
(LUE) whose treatment is very detailed. Correlation functions for two more matrix 
models - Jacobi Unitary Ensemble (JUE) and Cauchy Unitary Ensemble (CyUE) - 
are addressed in the Appendices C and D. 

Finally, in Section 5, we apply the integrable theory of CFCP to a comparative 
analysis of three alternative formulations of the replica method, with a special emphasis 
placed on the phenomenon of fermionic-bosonic factorisation of spectral correlation 
functions; some technical calculations involving functions of parabolic cylinder are 
collected in Appendix F. To make the paper self-sufficient, we have included Appendix 
E containing an overview of very basic facts on the six Painleve transcendents and a 
closely related differential equation belonging to the Chazy I class. 

The conclusions are presented in Section 6. 

2. Structure of the Theory 

The correlation function 



to be considered in this section can be viewed as a natural extension of its primary 
definition Eq. (1.1). Written in the eigenvalue representation (induced by the unitary 
rotation 5C = U^AIX such that U G U{n) and A — diag(Ai, • • • , A„)) it accommodates 
an n-dependent confinement potential * V^(A) and also allows for a generic eigenvalue 
integration domain f 



^ The properties of Hirota differential operators are reviewed in Appendix B. 
+ An alternative derivation can be found in Appendix A. 

* Matrix integrals with n-dependent weights are known to appear in the bosonic formulations of 
replica field theories, sec Osipov and Kanzieper (2007). This was the motivation behind our choice of 
the definition Eq. (2.1). 

t In applications to be considered in Section 5, the integration domain D will be set to D = [—1, +1] 
for (compact) fermionic replicas, and to D = [0, +00] for (noncompact) bosonic replicas. A more 
general setting Eq. (2.2) does not complicate the theory we present. 




(2.1) 



r 




(2.2) 
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The normalisation constant K„ is 

^ n 

^n= / T\dXje-''-^^'^-Al{\). (2.3) 

While for k G Z^., the correlation function n„|p(<j;K;) can readily be calculated by 
utilising the formalism due to Fyodorov and Strahov (2003), there seems to be no 
simple extension of their method to «; S M^. It is this latter domain that will be 
covered by our theory. 

Contrary to the existing approaches which represent the CFCP explicitly in a 
determinant form (akin to Eq. (1.12)), our formalism does not yield any closed 
expression for the correlation function n„|p(^;/«). Instead, it describes n„|p(^;K) 
implicitly in terms of a solution to a nonlinear (partial) differential equation which 
- along with an infinite set of nonlinear (partial) differential hierarchies satisfied by 
nn|p(^;K) - can be generated in a regular way starting with Eq. (2.1). Let us stress 
that a lack of explicitness is by no means a weak point of our theory: the representations 
emerging from it save the day when a replica limit is implemented in Eq. (1.9). 

Before plunging into the technicalities of the integrable theory of CFCP, we wish to 
outline its general structure. 

Deformation. — To determine the correlation function n„|p(<?; k.) nonperturbatively, we 
adopt the "dcform-and-study" approach, a standard string theory method of revealing 
hidden structures. Its main idea consists of "embedding" n„|p(^; k) into a more general 
theory of the r function 

r(«)(^,^.;t) = - / nK rn-.(^,'*;A,)e"(*^^^))-A^(A) (2.4) 

which posseses an infinite-dimensional parameter space (s;t) = (s; ti, t2, • " ") arising as 
the result of the (s; t)-deformation of the weight function 

r„(<^, «; A) = e-^"(^) n (fe - A)«» (2.5) 

Q = l 

appearing in the original definition Eq. (2.1). The parameter s is assumed to be an 
integer, s e Z, and v{t\ A) is defined as an infinite series 

oo 

t;(t;A) = ^^feA^ t = (ti, ta, • • •)• (2-6) 
fc=i 

Notice that a somewhat unusual (s G Z)-deformation of r„(<;, k;A) is needed to ac- 
count for the n-dependent confinement potential V^(A) in Eq. (2.1). 

Bilinear identity and integrable hierarchies. — Having embedded the correlation func- 
tion n„|p(^;/«) into a set of t functions Tn\<i, n]t), one studies the evolution of r 
functions in the extended parameter space (n, s, t) in order to identify nontrivial non- 
linear differential hierarchical relations between them. It turns out that an infinite set of 
hierarchically structured nonlinear differential equations in the variables t = {ti,t2,- ■ ■) 
can be encoded into a single bilinear identity 

I r(?)(t' - [z-^]) "^+^'"^.'^'/l+^""^^ 

I rf,e-(^-*'^-)rf+---^(t-[z-^]) ^-+^'^i;^+J"'^\ (2.7) 

Coo ^ 
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where the integration contour Goo encompasses the point z = oo. Here, a e M is 
a free parameter; the notation t ± [z~^] stands for the infinite set of parameters 
{tk ± /k}kez+', for brevity, the physical parameters <; and k were dropped from 
the arguments of r functions. 

To the best of our knowledge, this is the most general form of bilinear identity 
that has ever appeared in the literature for Hermitean matrix models: not only it 
accounts for the n-dependent probability measure ("confinement potential") but also 
it generates, in a unified way, a whole zoo of intcgrable hierarchies satisfied by the t 
function Eq. (2.4). The latter was made possible by the daedal introduction of the free 
parameter a in Eq. (2.7) prompted by the study by Tu, Shaw and Yen (1996). 

The bilinear identity generates various integrable hierarchies in the (n, s, t) space. 
The Kadomtsev-Petviashvili (KP) hierarchy, 

\ D,Du ri^^ it) o r(^) {t) = s^+.m) r^'\t) o r(^)(t), (2.8) 
and the Toda Lattice (TL) hierarchy, 

\ D^D, r(«) (t) o ri^) (t) = ([£»]) r^;^^ (*) ° rti^ (2-9) 

are central to our approach. In the above formulae, the vector D stands for D = 
(£>!, D2, • • • , -Dfe, • • •) whilst the fc-th component of the vector [D] equals k~^Dk. The 
operator symbol Dk f{t) o g{t) denotes the Hirota derivative * 

Dkf{t)og(t)^^f{t + x)g{t-x) . (2.10) 

The functions Sk{t) are the Schur polynomials (Macdonald 1998) defined by the 
expansion 

(00 \ 00 

E^j^' =E^'^^(*)' (2.11) 

see also Table 1. A complete list of emerging hierarchies will be presented in Section 
3.4. 

Projection. — The projection formula 

nn|p(«^;«) = ^T(«H«^,«;*) (2.12) 

J\ji s— 0, £—0 

makes it tempting to assume that nonlinear integrable hierarchies satisfied by r func- 
tions in the (n, s, f )-spacc should induce similar, hierarchically structured, nonlinear 
differential equations for the correlation function n„|p(ii;K;). To identify them, one 
has to seek an additional block of the theory that would make a link between the 
partial {tk}kei+ derivatives of r fimctions taken at t = and the partial derivatives 
of n„|p(<j;K;) over the physical parameters {Ta}aez+- The study by Adler, Shiota and 
van Moerbeke (1995) suggests that the missing block is the Virasoro constraints for r 
functions. 

Virasoro constraints. — The Virasoro constraints reflect the invariance of r functions 
[Eq. (2.4)] under a change of integration variables. In the context of CFCP, it is useful 
to demand the invariance under an infinite set of transformations * 

o 

A.^^. + e^f V(M,)n(M,-4), q>-l, (2.13) 

k=l 

* The properties of Hirota differential operators are briefly reviewed in Appendix B; see also the book 
by Hirota (2004). 

* The specific choice Eq. (2.13) will be advocated in Section 3.6. 
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Table 1. Explicit formulae for the lowest-order Schur polynomials s^(t) defined 
by the relation exp (E^^i tjxA = T.'^o 



e 


i\s 


i{t) 









1 








1 


ti 








2 


tl- 


f 2*2 






3 


tl- 


f 6tlt2 +6*3 






4 


ti- 


f 24*1*3 + 12*?*2 - 


h 12*1 + 


24*4 


5 


4- 


f 20*f *2 + 60*f *3 - 


h 60*1*1 


+ 120*1*4 + 120*2*3 + 120*6 



The Schur polynomials admit the representation (Macdonald 1998) 

where the summation runs over all partitions A = (£j ^ , ■ ■ ■ , f-g^ ) of the size | A| = i. 
The notation A = ((.'^^ , ■ ■ ■ , ig^ ), known as the frequency representation of the 
partition A of the size |A| = £, implies that the part ij appears aj times so that 
i = "^j—i^jf^j, where g is the number of inequivalent parts of the partition. 
Another way to compute S({t) is based on the recursion equation 
1 ^ 

3 = 1 

supplemented by the condition so(t) = 1- 



labeled by integers q. Here, e > 0, the vector c' is c' = {ci, • • • , C2r} \ {±00, Zq} with 
Zq denoting a set of zeros of /(A), and g = dini(c'). The function /(A) is, in turn, 
related to the confinement potential V^„_s(A) through the paramcterisation 

%^ = 77^' ^w = f:M^ m = f:a,x'' (2.14) 

in which both g{X) and /(A) depend on n—s as do the coefficients bk and au in the above 
expansions. The transformation Eq. (2.13) induces the Virasoro-Uke constraints jt 



where the differential operators 

e 



C^{t)+lf\<;;t)\ rW(^;t) = :Bj^(0 rW(«^; *), (2-15) 
tial operators 

(*) = E E se-k{-p,ic')) ( '^e^,+k+i{t) - h- (2.16) 

e=Q fe=o ^ q+k+e+i / 

and 

00 g q+k+e / p \ „ 

Lf\t) = Y^a,Y.s,^,{-p^ic')) Y: E^"C L, ^ ^ (2.17) 

i=0 fe=0 m=0 \a=l J Olq+k+t-m 

act in the t-space whilst the differential operator 



= E ( n(^« - 4) ) /(ca)cr^^ (2.18) 

a=l \fe=l 



9Ca 



tt The very notation £^ suggests that this operator originates from the confinement-potential-part 

g-V,! in Eqs. (2.5) and (2.4). On the contrary, the operator is due to the determinant-like product 

na('>« ~ A)""* in Eq. (2.5). Indeed, setting = nullifies the operator C'^"^. See Section 3.6 for a 
detailed derivation. 
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acts in the space of physical parameters {<?a}aez+- The notation Sfc(— p^(c')) stands 
for the Schur polynomial and Pg{c') is an infinite dimensional vector 



p,(c') = (tr,(c'), \tr,{c')\ ^tr,(c')^ • • •) (2.19) 



with 



tr,(c')'= = l^(c;)^ (2.20) 
Notice that the operator jt^{t) is expressed in terms of the Virasoro operators ff 



obeying the Virasoro algebra 

[Lp,Lg] = {p-q)Lp+g, p,q>-l. (2.22) 



Projection (continued). Equations (2.12) and (2.15) suggest that there exists an 
infinite set of equations which express various combinations of the derivatives 



d 



dt, 



,2 



s=a.t=0 



in terms of !]3^(^) n„|j,(^; k.). This observation makes it tempting to project the 
hierarchical relations Eqs. (2.8) and (2.9) onto the hyperplane (s = 0,t = 0) in an 

attempt to generate their analogues in the space of physical parameters. In particular, 
such a projection of the first equation of the KP hierarchy, 

+ 4 - "si) ^""'<" + " (5 '"^ '^"'")' = 

is expected § to bring a closed nonlinear differential equation for the correlation function 
nn|p(<>; 1^) of characteristic polynomials. It is this equation which, being supplemented 
by appropriate boundary conditions, provides an exact, nonperturbative description 
of the averages of characteristic polynomials. Similarly, projections of other equations 
from the hierarchies Eqs. (2.8) and (2.9) will reveal additional nontrivial nonlinear 
differential relations that would involve not only n„|p(<j; k) but its "neighbours" 
n„-|-q|p(i;; k), as explained in Section 4. 

Having exhibited the general structure of the theory, let us turn to the detailed 
exposition of its main ingredients. 



3. From Chgiracteristic Polynomials to r Functions 

3.1. The T function, symmetry and associated scalar product 

Intcgrability derives from the symmetry. In the context of t fimctions Eq. (2.4), the 
symmetry is encoded into A^(A), the squared Vandermonde determinant, as it appears 
in the integrand below ||: 

^n\t) = ^ / n id\,Tn-s{\i)e<'-'^A-^l{X). (3.1) 

ft For q = —1, the second sum in Eq. (2.21) is interpreted as zero. 

§ Whether or not the projected Virasoro constraints and the hierarchical equations always form a 
closed system is a separate question that lies beyond the scope of the present paper. 
II For the saJce of brevity, the physical parameters <i and n were dropped from the arguments of r„ 
and r„_3. 
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In the random matrix theory language, the r function Eq. (3.1) is said to posses the 
13 = 2 symmetry. Using the identity 

A„(A) = det[A^^-i]i<,- fe<„ = det[Pfe_i(A,)]i<,-,fc<„ (3.2) 

with Pfe(A) being an arbitrary set of monic polynomials and the integration formula 
(Andreief 1883, de Bruijn 1955) 



(A)Vfc(A) 



the T function Eq. (3.1) can be written as the determinant 

r(«)(*)=det Ur)(<)] 
of the matrix of moments 

M$r^(*) = (^^.l^'fe)r.e«= / rfAr„(A)e^(*^^)p,(A)P,(A) 
Both the determinant representation and the scalar product 

{f\9)^= [ dXw{X)mg{X) 
are dictated by the j3 = 2 symmetry § of the r function. 



(3.3) 

(3.4) 
(3.5) 
(3.6) 



3.2. Bilinear identity in integral form 

In this subsection, the bilinear identity Eq. (2.7) will be proven. 

The T function and orthogonal polynomials. — The representation Eq. (3.4) reveals a 
special role played by the monic polynomials Pl"^\t; A) orthogonal on D with respect 
to the mecisure ry7^(A) e^^^'^-^dA. Indeed, the orthogonality relation 

{Pk\Pj)r^,. = f rfAr„(A)e''(*^^)pf )(t; A)P,("')(t; A) = ht\t)6jk, (3.7) 

shows that the choice Pj{X) i-> Pj"~*''(t;A) diagonalises the matrix of moments in 
Eq. (3.4) resulting in the fairly compact representation 

n-l 

^^^^W^H'^r^^W- (3-8) 

§ The T function Eq. (3.1) is a particulcir case of a more general r function 

1 /■ " 

In accordance with the Dyson "tlircc-fold" way (Dyson 1962), the symmetry parameter /3 may also 
take the values /3 = 1 and (3 = 4. For tlicsc cases, the t function Eq. (3.7) admits the Pfaffian rather 
than determinant representation (Adlcr and van Moerbeke 2001); 

.«(*;,)= pf Hr^' 

where the matrix of moments /i^.^'(t;/3) = {Pj\Pk)^^ gv is defined through the scalar product 



{f\9) 



(/3) 



[ dXdX'w(X) /(A) sgn (A' - X)w(X')g(X'), 

71)2 

f dXw{X) [/(A)g'(A)-ff(A)/'(A)] , 



1; 

: 4. 
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Remarkably, the monic orthogonal polynomials P^'"^(t; A), that were introduced as 
a most economic tool for the calculation of Tn^ , can themselves be expressed in terms 
of T functions: 



Prit;X) = X^ (3.9) 

Here, the notation t — [\~^] stands for an infinite-dimensional vector with the 
components 

t ± [\-'] = (^ti ± ^, i2 ± • • • , tfc ± • • •) . (3.10) 

The statement Eq. (3.9) readily follows from the definitions Eqs. (3.1) and (2.6), the 
formal relation 

g.(t±[A-M; A,) ^ g.(t; A,) -^y\ (3.11) 
and the Heine formula (Heine 1878, Szego 1939) 

Pt\t;X)= [ [ n(dA,(A-A,)r„(A,)e''(*^^^)).A2(A). (3.12) 



The T function and Cauchy transform of orthogonal polynomials. — As will be seen 
later, the Cauchy transform of orthogonal polynomials is an important ingredient of 
our proof of the bilinear identity. Viewed as the scalar product, 

Qi"'\t;z) = (Pt\t;x) = / dAr„(A)e''(*^^) (3.13) 

it can also be expressed in terms of r function: 

(fc+l-m) /, r-,-ll\ 

Qi-\t;.)=z->'-^ ^^-^ J^ + f (3.14) 

Tk (*) 

To prove Eq. (3.14), we substitute Eq. (3.12) into Eq. (3.13) to derive: 



{z — Xk+i ) j-J^ Xk+1 — A_ 
Owing to the identity 



" 1 " / 1 " 1 \ 



taken at n = A; + 1, the factor 
1 ^ 



{z - Xk+i) j-J^ Xk+i - Xj 
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in the integrand of Eq. (3.15) can be symmetrised 

1 A 1 ^, 1 ^/^ 1 'f^ 1 



a=l \ j=l, i^a 

fe+1 



{z — A/s+i) Afe+i — Aj fc 

fc + 1 -j-J- ^ - A,- 

to yield the representation 



In view of Eq. (3.11), this is seen to coincide with 

^-fe-i I- fe+i 

(fc + 1)!,, 

Comparison with the definition Eq. (3.1) completes the proof of Eq. (3.14). 



Proof of the bilinear identity. — Now we are ready to prove the bilinear identity 

U+l+s-m),. r -11-, 

,,e(°-^)-(^-^'^-)ri-)(t'-[z-^]) "^+^ + 

i dz e"''^*-^'--) rf (t - [z-^]) '^i^^t^^"'^^ , (3.18) 

Coo 

where the integration contour Coo encompasses the point ^ = oo, and a € M is a free 
parameter. 

We start with the needlessly fancy identity 
/ dAr„(A) e''(*;^)e(«-i)''(*-*''^)p,(")(t; A)Pi")(*'; A) 

= [ rfAr„(A)e''(*''^)e<"'(*-*''^)p^(")(t;A)Pi")(t';A) (3.19) 

whose structure is prompted by the scalar product Eq. (3.7) and which trivially holds 
due to a linearity of the i-deformation 

v{t; A) + (a - 1) v{t - t'; A) = v{t'; A) +av{t- t'; A), (3.20) 

see Eq. (2.6). 

The formulae relating the orthogonal polynomials and their Cauchy transforms to 
T functions [Eqs. (3.9) and (3.14)] make it possible to express both sides of Eq. (3.19) 
in terms of r functions with shifted arguments. 

(i) Due to the Cauchy integral representation 

g(a-lMt-t';A)pH(f';A) = -L / ^^e(a-lMt-t';.):^|!lM, (3.21) 

the l.h.s. of Eq. (3.19) can be transformed as follows: 

l.h.s. = — / dze(''-iM*-*';^)Pi")(t';z) / dAr„(A)e-(*;^):^£i^ 
2Tn Jq^ z-X 



Q[''\t;z) [Eq. (3.13)] 



dze^"-^^''^*-*'''^ Pj;^\t';z)Q'f-\t;z). (3.22) 
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Taking into account Eqs. (3.9) and (3.14), this is further reduced to 
l.h.s.= ' 



2'rirf-"^(t)Ti"-")(t') 

x/ rf.e("-^)"(*-*'^-) rt-Ht' - [^-^]) "^+"^"ltt^""^^ . (3.23) 



(ii) To transform the r.h.s. of Eq. (3.19), we make use of the Cauchy theorem in 
the form 

^av(t-t';X)pin)^^,^^^J_ f ^av{t-t' ;z) ^l^}^ ^ (334) 

to get: 

r.h.s. = -^(/' dze''<*-'''^^P^^\t-,z) [ dAr„(A)e''<*''^) ^'^"^^*;^^ 

QW(t';z) [Eq. (3.13)] 

= dze'^^('-''-'^^Pt\t;z)Q^;^\t';z). (3.25) 

Taking into account Eqs. (3.9) and (3.14), this is further reduced to 
r.h.s. = 



27rirf-")(f)Ti'"-")(t') 



(3.26) 



The bihncar identity Eq. (3.18) follows from Eqs. (3.23) and (3.26) after setting 
n = m — s. End of proof. 

3.3. Bilinear identity in Hirota form 

The bilinear identity Eq. (3.18) can alternatively be written in the Hirota form: 

oo 

e^(-^)^s,((2a-l-/3)x)s,+,+i([i5]) rj^'Ht) o r^l+'^\t) 

k=0 

oo 

^^-Pi.-m J2 5,((2a-l + /3)x)a._,_i([i^])r(^+^+^)(t)or(r/)(t) (3.27) 

k=q+l 

where /3 = ±1 (not to be confused with the Dyson symmetry index!), p > 1 
and q > —1. Let us remind that the vector D appearing in the scalar product 
(a; • D) = ^j.XkDk stands for D = {Di, D2, ■ ■ ■ , Dk, ■ ' ')] the fc-th component of 
the vector [D] equals k^^Dk (compare with Eq. (3.10)). The generic Hirota differ- 
ential operator T(£>) f{t) o g[t) is defined in Appendix B. Also, Sk{x) are the Schur 
polynomials defined in Eq. (2.11). 

To prove Eq. (3.27), we proceed in two steps. 

(i) First, we set the vectors t and t' in Eq. (3.18) to be 

{t,t')^ {t + x,t-x). (3.28) 

The parameterisation Eq. (3.28) allows us to rewrite the (t, t') dependent part of the 
integrand in the l.h.s. of Eq. (3.18) 
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as follows: 

exp [v{2ia - l)x; z)] {t-x- [z-'])4lX'+'-"'^ it + x + [z'^]) 
= exp [w(2(a — l)x; z)] 

X exp [{X . d^) + i[z-'] . d^)] T^if+'-'Ht + ri^Ht - (3-29) 

Here, we have used the linearity of the f-deformation, a v{t; z) = v{at; z). Further, we 
spot the identity ([z"^] ■ d{) = v {[d^]; z''^) to reduce Eq. (3.29) to % 

exp [t;(2(a — l)x; 2;)] 

X exp [{x ■ d^)] exp [v ([9^; z-')] t\%'+'-"'^ [t + i) r^^ [t - 

The latter can be rewritten in terms of Hirota differential operators [see Eq. (B.2)] with 
the result being 

exp [v{2{a - l)x; z)] exp [{x ■ D)] exp [v {[D]- z-^)] Tf+^^+'-^\t) o Ti«)(t). (3.30) 
By the same token, the (t, t') dependent part of the integrand in the r.h.s. of Eq. (3.18), 

can be reduced to 

exp[^;(2a^r;z)]exph(a;.D)] exp ([£»]; z-^)] riH'^ {t) o Tf+'-^\t). (3.31) 
Thus, wc end up with the alternative representation for the bilinear identity Eq. (3.18): 

e^^-""^/ 7I^exp[^;(2(a-l)a=;.)] e^v[vm,z-')]rf,\'^'--\t) or^^){t) 

= e"^^-''^ / 3^+1 exp[.(2ax;.)] exp [. ([D]; ."!)] r^^Xl\t) o rf^^-\t). 

(3.32) 

(ii) Second, to facilitate the integration in Eq. (3.32), we rewrite the integrands 
therein in the form of Laurent series in z by employing the identity 



= expK^t,z-'- =J2sk{t)z\ (3.33) 
/ fe=o 

Now, the integrals in Eq. (3.32) are easily performed to yield 

oo 

e(-^) ^ s, (2(a - l)x) s^+m-i {[D]) ri^f ^^"'"^ {t) o r(j) (t) 

fc— max(0,^ — m) 
oo 

= J2 Sk{2ax)sk+i-m{[D])T^^+l\t)oT^/+'-"'\t). (3.34) 

fc=max(0,m— ^) 

It remains to verify that Eq. (3.34) is equivalent to the announced result Eq. (3.27). 

To this end we distinguish between two different cases, (i) If £ < m, wc set £ = p — 1, 
m = p + q and s i-> s + g in Eq. (3.34) to find out that it reduces to Eq. (3.27) taken 
at /3 = +1; (ii) If ^ > m, we set ^ = p + m = p — 1 and ,s s — 1 in Eq. (3.34) to 
find out that it reduces to Eq. (3.27) taken at /3 = — 1. This ends the proof. 

For an alternative derivation of Eq. (3.27) the reader is referred to Appendix A. 

t Here, 

, /did Id 
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3.4- Zoo of integrable hierarchies 

The bilinear identity, in either form, encodes an infinite set of hierarchically structured 
nonlinear differential equations in the variables t. Two of these hierarchies — the KP 
and the TL hierarchies - were mentioned in Section 2. Below, we provide a complete 
list of integrable hierarchies associated with the r function Eq. (3.1). 

To identify them, we expand tlic^ bilinear identity in Hirota form [Eq. (3.27)] around 
a; = and a = 0, keeping only linear in x terms. Since so{t) = 1 and 

Sk{t) =tk + 0{t''), /e = l, 2, ... (3.35) 

we obtain: 

(l-<5,,_i)s,+i([r>])r(«)(t)or(f/(t) 

oo 

+ J2^k [(2a-l-/3)sfe+,+i([r>])+/3A(s,+i([-D])+^,,-i)] t(«)(*) o^Z'^*) 

k=l 

OO 

-(2a-l + /3) Yl a;feSfe_,_i([i?])T(;+jV^(t)or^!.-i^)(t) + C)(x2) = 0. (3.36) 

fc— max(l,g+l) 

As soon as Eq. (3.27) holds for arbitrary a and x, Eq. (3.36) generates four identities. 

(i) The first identity 

im r(«)(t) o r(;+/(t) = (3.37) 
holds for g' > 1 and k = 0, 1, . . . , q. 

(ii) The second identity 

[(1 + /3) im - l3DkS,+, i[D]) ] r(^)(t) o r^l+'^\t) = (3.38) 

holds for g > 1 and A; = 1, 2, . . . , g. 

(iii) The third identity 

[(1 + (3) Sk+,+i m) - PDu{s,+, m) + 5,,_i)] r^^)(t) o r^;+/)(t) 

= (1 - /3) ([£>]) t(;+^+^) (t) o rl^S,'\t) (3.39) 

holds for g > — 1 and k > max(l, q+1). 

(iv) The last, fourth identity 

m) 4'\t) o = m) rizt^'\t) o Ti:_-,'\t) (3.40) 

holds for g > and k>q + l. 

Equations (3.37) - (3.40) can further be classified to yield the following bilinear hi- 
erarchies: 



• Toda Lattice (TL) hierarchy: 

\d,D, r(^) {t) o r(«) (t) = sk-, m) r(;V^ W ° i-i'^ (*) (3-41) 
with A; > 1. 
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q-modified Toda Lattice hierarchy: 
1 

2' 

with g > and k > q + 1. (For g = 0, it reduces to the above Toda Lattice 
hierarchy.) 



Kadomtsev-Petviashvili (KP) hierarchy: 

^D^Dk - Sk+i {[D]) 
with t) fc > 3. 



°^i''^(*)=0 (3.43) 



q-modified Kadomtsev-Petviashvili hierarchy: 
1 
2 

with q>0 and fc > g + 1. (For g = 0, it reduces to the above KP hierarchy.) 



{[D]) TW(t) o r(;+/^(t) = ([r>]) rW(t) o r^;+;)(t) (3.44) 



• Left q-modified Kadomtsev-Petviashvili hierarchy: 

D,s,+, im r(^)(t) o r(f/)(t) = (3.45) 
with g > 1 and 1 < < g. 

• Right q-modified Kadomtsev-Petviashvili hierarchy: 

([£>]) TW(t) o r^;t'^(t) = (3.46) 
with g > 1 and < < g. 

• (—l)-modified Kadomtsev-Petviashvili hierarchy: 

h ([£)])- 7^,] r(^)(t) or^!_-'\t)=0 (3.47) 
with k>2. 

Notice, that the modified hierarchies will play no role in further consideration. 

3.5. KP and Toda Lattice hierarchies 

As was pointed out in Section 2, the KP and Toda Lattice hierarchies are of primary 
importance for our formalism. In this subsection, we explicitly present a few first mem- 
bers of these hierarchies. 

KP hierarchy. — Due to the properties of Hirota symbol reviewed in Appendix B, the 
first nontrivial equation of the KP hierarchy corresponds to fc = 3 in Eq. (3.43). Con- 
sulting Table 1 and having in mind that [D]k = k~^Dk, we derive the first two mem- 
bers, KPi and KP2, of the KP hierarchy in Hirota form 

KPi : {Df - 4D,Ds + 3DI) r^«) (t) o t^«) (t) = 0, (3.48) 

KP2 : {DID2 + 2D2D3 - 3DiDi) r^«) (t) o r^«) (t) = 0. (3.49) 

In deriving Eqs. (3.48) and (3.49), we have used the Property 2a from Appendix B. 

tl Both fe = 1 and k = 2 bring trivial statements, see Appendix B. 



3 From Characteristic Polynomials to r Functions 



19 



Making use of the Property 2b from Appendix B, the two equations can be written 
expUcitly: 



Vaif5i2 dtidt4 at25t3, 

+6(|log.(^)w)(^log.(^)w)=0. (3.51) 



Only KPi will further be used. 



TU : ^Dl r(«)(t) o r(«)(t) = r(:^^)(t) o r^;_-'\t), (3.52) 



Torfa Lattice hierarchy. — The first nontrivial equations of the Toda Lattice hierarchy 
can be derived along the same lines from Eq. (3.41). 

1 

2^ 

TU : Id,D, rW(t) o r(«)(t) = D, 4lf\t) o 4'_-,'\t). (3.53) 
Explicitly, one has: 

TLi : '{t) — I I = r^+i '(t)T^-i '{t), (3.54) 

dT, dlT 

= r(l-^(t)^^£i^-ri^^^^ (3.55) 

Higher order members of the KP and Toda Lattice hierarchies can readily be 
generated from Eqs. (3.43) and (3.41), respectively. 



3.6. Virasoro constraints 

Virasoro constraints satisfied by the t function Eq. (3.56) below is yet another 
important ingredient of the "deform-and-study" approach to the correlation functions 
of characteristic polynomials n„|p(<j; k). In accordance with the discussion in Section 2, 
Virasoro constraints are needed to translate nonlinear integrable hierarchies Eqs. (3.41) 
- (3.47), satisfied by the r function 

TW(<^,K;t) = l / n (dA,e-^— (^^) n(^„-A,)««e''(*^^^)) •A2(A), (3.56) 

into nonlinear, hierarchically structured differential equations for the correlation 
function 

n„|p(^; n) = ^ J^^ fl ^dXj e-^"(^^) (^„ - A,)- j • A^(A) (3.57) 

obtained from Eq. (3.56) by setting t = and s = 0. The normalisation constant J^n 
is defined in Eq. (2.3). 

The Virasoro constraints reflect invariance of the r function Eq. (3.56) under the 
change of integration variables 

Xj iij + enf^^Riiij), q > -1, (3.58) 
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labeled by the integer q; here e > is an infinitesimally small parameter, and R{fJ,) 
is a suitable benign function (e.g., a polynomial). The function /(A) is related to the 
confinement potential Vn-sW through the parameterisation (Adler, Shiota and van 
Moerbeke 1995) 



dVn- 



dX 



5(A) 
/(A) 



(3.59) 



k=0 fe=0 

in which both g{X) and /(A) depend on n — s as do the coefficients bk and Ofc in the 
above expansions. We also assume that 

lim /(A) A'' e-^"-=(^) =0, k>Q. (3.60) 

A->-±cx> 

To derive the Virasoro constraints announced in Eqs. (2.15)-(2.18), we transform 
the integration variables in Eq. (3.56) as specified in Eq. (3.58) and further expand 
Eq. (3.56) in e. Invariance of the integral under this transformation implies that the 
linear in e terms must vanish: 

/ (..)f,E^^''^"^"'^"''^'^^ " 



Hi - iij 



{q + l)lJ,l+v'{t;iJ,e)lJ.i 



9+1 



(dim(c') \ 



Here, c' = {ci, • • • , C2r} \ {±00}, 

n 

fe=i 
and 

n 



c«=l 



(3.61) 



(3.62) 



(3.63) 



fc=i 



In the above formulae, we reinstated /3 > 0; it will be set to ^ = 2 when needed. 

The choice of R{^j) is dictated by problem in question and, hence, is not unique. 
If one is interested in studying matrix integrals as functions of the parameters 
{ci, • • • , C2r} defining the integration domain D, the suitable choice of R{ij) is 

R{^i) = /(/x). (3.64) 

In this case, the differential operator (Adler, Shiota and van Moerbeke 1995) 

Y^^iOc^r'^ (3.65) 

becomes an essential part of the Virasoro constraints. In the context of characteristic 
polynomials, the integration domain D is normally fixed whilst the physical parameters 
{<Jct} are allowed to vary. This prompts the choice 



R{l^) = /(a*) fl (/^ - 4), Q = dim(c') 



(3.66) 



fe=i 
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that nullifies the differential operator Eq. (3.65). Equivalently, 

o 

= /(m) (3-67) 



fc=0 



Here, the notation Sk{—Pg{c')) stands for the Schur polynomial and Pg{c') is an infinite 
dimensional vector 

Pg{c') = (tr,(c'), ^tr,(c')^ • • • , ^tr,(c')^ • • •) (3.68) 

with 

tr,(c')'==X:(c;)^ (3.69) 

Rem,ark. Equation (3.66) assumes that none of c^'s are zeros of /(/u). If this is not the 
case, the set c' must be redefined: 

c'^c'\{Zo}, (3.70) 

where Z.o is comprised of zeros of /(/i). 

Substituting Eqs. (3.67), (3.59) and (2.6) into Eq. (3.61), we derive: 

oo / „ q+k+i 

E«^(f E tr„(M^)tr„()u''+'=+'-^) 



+ {i + k + q+l) tr„(/x«+'=+^) + Ei*itr„(/x«+'=+'+^'^ 



p q+k+i p „ N 

+E«" E ctr„(iu«+'=+^-™)-E^r'=+^+^^ 

a=l m=0 a=l y 

e 



(d/x) E^^-'fe(-P.(^'))E^^tr„(/x«+'=+^+^)4«Ht;/x,0. (3.71) 

fe=0 i=0 

The ^-dependent part in Eq. (3.71), 
/ {dn)^Sg-ki-Pe{c'))^a, 

(p g+fe+i p o \ 

E«« E Ctr„(/x''+'=+^-'")-E^«^'^'^V ^n\t;i^,<;), (3.72) 

originates from the term 

-/ id^.)(j2^^r'R{^^e)j27^]^^nHt■,^^,<;) (3.73) 

in Eq. (3.61). Indeed, substituting Eqs. (3.67) and (3.59) into Eq. (3.73), the latter 
reduces to 

/Q oo / p n q+k+i+1 \ 

(dn) E«.-fc(-p.(c'))E«' E^«E _ ii'Ht;n,<^)- (3.74) 

The double sum in parentheses can conveniently be divided into two pieces, 

P " O + fe+i+l g + fc + j+l 
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and 



P " q+k+i+l 

Due to the identities 

n q+k+i+l _ q+k+i+l Q+k+i 

yffi \. ^ y ^^^^^^^q+k+i-m^ (377) 

and 

E (*; = - ^4'^ (*; (3.78) 

we conclude that Eq. (3.73) reduces to the sought Eq. (3.72). We found it more 
convenient to rewrite the d/d<;a-teTm in Eq. (3.72) in a more compact way, 

•^D" fe=0 i=0 Va=l "^^"Z 

= ^^-(0 / {d^^)I(c\t■,^l,<;) (3.79) 

with the differential operator 'B^(^) being 

= E ( n('^" -4)) /(?a)cr^^. (3.80) 

Equation (3.79) follows from the expansions Eqs. (3.59), (3.66) and (3.67). 

To complete the derivation of Virasoro constraints, wc further notice that terms 
tr„(ju-') in Eq. (3.71) can be generated by differentiating In ^ over tj. Since tr„(/LX°) = n, 
the derivative d/dto should formally be understood as d/dto = n. This observation 
yields Virasoro constraints for the r function 

(^, t) = 1 / {d,x) (t ; /X, ^) (3.81) 

in the form {q > —1) 

£j'(f) +£f '(^;t)] r(^)(^;t) = ®^(^) r^^) (^; t). (3.82) 
Here, the differential operators 

^J'W = E E ^.--^(-P.CcO) («^>cf^),+,(t) - h,^^——^^ (3.83) 



£=0 fe=0 

and 

oo Q q+k+e / p 



d 



£f*(t) = ^a,^s,_fe(-p,(c')) E«"CUt (3.84) 

e=0 fe=0 m=0 Vc«=l / (^Q+k+i-m 

act in the t-space whilst the differential operator ^ J'(^) acts in the space of physical 
parameters {<;a}a<£Z+- Notice that the operator 'C^(t) is expressed in terms of the 
Virasoro operators 

^f«=g^'.^4i:^^0-f)(»-^|;. (3-) 

obeying the Virasoro algebra 

[^i^\^f^] = {p-Q)^^J2q, P,9>-1- (3.86) 
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The Virasoro constraints derived in this section stay vaUd for arbitrary ^ > 0; for 
/3 = 2, they are reduced to Eqs. (2.15) - (2.22) announced in Sec. 2. 

This concludes the derivation of three main ingredients of integrable theory of aver- 
age characteristic polynomials - the bilinear identity, integrable hierarchies emanating 
from it, and the Virasoro constraints. 



4. Prom T Functions to Characteristic Polynomials 

The general calculational scheme formulated in Section 2 and detailed in Section 3 
applies to a variety of random matrix ensembles. In this Section we deal with CFCP 
for the Gaussian Unitary Ensemble (GUE) and Laguerre Unitary Ensemble (LUE). A 
detailed treatment of the GUE case is needed to lay the basis for further comparative 
analysis of the throe variations of the replica approach that will be presented in Section 

5. The study of the LUE relevant to the QCD physics (Verbaarschot 2010) is included 
for didactic purposes. A sketchy exposition of the theory for Jacobi Unitary Ensemble 
(JUE) and Cauchy Unitary Ensemble (CyUE) appearing in the context of universal 
aspects of quantum transport through chaotic cavities (Beenakker 1997) can be found 
in Appendices C and D. 

4-1- Gaussian Unitary Ensemble (GUE) 

The correlation function of characteristic polynomials in GUE is defined by the n-fold 
integral 

n^ip(^; '^^=^ [ n r^^- n - ^^•)''" ) • (4.1) 



where 



11 n 

/ n ^"'0 • = w22-"("-i)/2 n + 1) (4-2) 

is the normalisation constant. The associated r function equals 

r^i^,K;t)^^ f fl(dX,e-'^+"^'-'^^ll{,^-X,)A-Al{X), (4.3) 

see Section 2. (In the above definitions, the superscript G stands for GUE but it will 
further be omitted when notational confusion is unlikely to arise.) 

4.1.1. Virasoro constraints 

In the notation of Section 3, the definition Eq. (4.1) implies that 

/(A) = l ^ ak = dk,o, (4.4) 

g{X)=2X ^ bk = 26k,i, (4.5) 

D = M ^ dim(c') = 0. (4.6) 

This brings the Virasoro constraints Eqs. (2.15) - (2.21) for the r function Eq. (4.3): 



r„(^,/t;t) = 23,r„(^,/t;t), (4.7) 
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where 

-I 



a=l 

the short-hand notation '&m{<i,K) is defined as 

p 

^„(^,k) = ^k„C, (4.9) 

a=l 

SO that 

p 

i?o(<», k) = trp K = = K. (4-10) 

Q = l 

Also, £5(4) is the Virasoro operator given by Eq. (2.21). Notice that the r function 
in Eq. (4.7) does not bear the superscript (s) since the GUE confinement potential 
V{\) = docs not depend on n. 

In what follows, we need the three lowest Virasoro constraints labeled by q = — 1, 
q = and q = +1. Written for the logarithm of r function, they read: 

00 Q d \ 

"^■^^^dt- '^df I = ^-1 logr„(^,K,;t), (4.11) 

X]i*J a7: - 2^ 1 logT„(^,K;t)-|-n(n + /t) =3o logT„(<j,K;t), (4.12) 



7=1 



2 



+ni?i(^,/«) = Si logr„(^,K;t). (4.13) 

Toda Lattice hierarchy 

Projection of the Toda Lattice hierarchy Eq. (2.9) for the t-dependent r function 
Eq. (4.3) onto the hyperplanc t = generates the Toda Lattice hierarchy for the 
correlation function n^|p(^; k) [Eq. (4.1)] of the GUE characteristic polynomials, 

Below, the first [Eq. (3.54)] and second [Eq. (3.55)] equation of the TL hierarchy will 
be considered: 

^^,logr„(t) ^^^^^ 

OT1OT2 T^(t) Oti \Tn-l{t) J 

The equivalence of Eqs. (4.15) and (4.16) to Eqs. (3.54) and (3.55) is easily established. 



Ti-i • "572 log^"(*) = 1277\ ' (4.15) 



(i) To derive the first equation of the TL hierarchy for n^|p(<j; k) from Eqs. (4.14) 
and (4.15), we have to determine 



& 
dtl 



2 

2 l0gT„(t) 



t=0 
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with the help of Virasoro constraints. This is achieved in two steps. First, we 
differentiate Eq. (4.11) over ti and set t = afterwards, to derive: 



2^1ogT„(t) 



n-'B-i — logT„(t) 



t=o 



Second, we set t = in Eq. (4.11) to identify the relation 



2Aiogr„(t) 
at I 



logr„(0). 



Combining Eqs. (4.17) and (4.18), we conclude that 

4^ logT„(t) = 2n + -Bl^ logT„(0). 

Finally, substituting this result back to Eq. (4.15), and taking into account Eqs. (4.14) 
and (4.2), we end up with the first TL equation 



(4.17) 
(4.18) 
(4.19) 



— G 

TLi : 



logn„|p(^;K) = 2n 



nn+i|p(^; *«) nn-ib(<?; «) 



- 1 



(4.20) 



written in the space of physical parameters ^. 



(ii) The second equation of the TL hierarchy for n„|p(^;/«) can be derived along the 
same lines. Equation (4.16) suggests that, in addition to the derivative d/dtilogTn at 
t = given by Eq. (4.18), one needs to know the mixed derivative 



52 



dtidt2 



logr„(t) 



t=o 



It can be calculated by combining Eq. (4.11) differentiated over t2 with Eqs. (4.12) and 
(4.18). The result reads: 



TLn 



(1-Bo):B_i logn„|p(.j;K) 

_ n„_|_i|p(^; k) n„_i|p(^; k) 



:B_ilog 



n 



n+l 



ip(<?;«^) 



n, 



n—l\p 



(4.21) 



Higher order equations of the TL hierarchy for the correlation functions Iln\p can be 
derived in a similar fashion. 

~G ~G 

Remark. — For p = 1, the equations TL^ and TLg become particularly simple: 

rpfG /n„+i((r;K)n„_i(c;K) \ , - 

TL, : ^logn„(,;.)=2n(^ ^ 1 j , (4.22) 

/ 9\ 9 ^ u ( \ n„+i(s;/,:)n„_Us:/v) 8 , f 

(i-.^j^iogn„(.;«) = n — — a^i°Hn;::T(^ 

(4.23) 

4.1.3. KP hierarchy and Painleve IV equation 

The technology used in the previous subsection can equally be employed to project 
the KP hirerachy Eq. (2.8) onto the hyperplane t = 0. Below, only the first KP 
equation 
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will be treated. Notice that no superscript (s) appears in Eq. (4.24) as the GUE 
confinement potential docs not depend on the matrix size n. Proceeding along the lines 
of the previous subsection, we make use of the three Virasoro constraints Eqs. (4.11) - 
(4.12) to derive: 

16^ logT„(t)|^^^ = ®ii logT„(0), (4.25) 
= 2n(n + K)-(2-Bo):Bo logT„(0), (4.26) 



log r„(t) 



t=0 



t=o 



= n (3n + 2k) 



1 



- ( ®o - ^1^-1 - § (2n + «) ^-1 ) logT„(0). 



(4.27) 



Substitution of Eqs. (4.25) - (4.27) and (4.17) into Eq. (4.24) generates a closed 
nonlinear differential equation for logn„|p(<j; k) in the form 



'h\ + 8(n - K)'h\ - 4(2«o - 3:Bg + ^M-x) log n„,p(^; k) 



+6 'Biilogn„|p(<;;K) = 8nK. 



(4.28) 



Notice that appearance of the single parameter k instead of the entire set k = 
{ki, . . . , Kp) in Eq. (4.28) indicates that correlation functions n„|p(^;K) with differ- 
ent K but with identical traces tr^ k satisfy the very same equation. It is the boundary 
conditions | that pick up the right solution for the given set k = (ki, . . . , Kp). 



Remark. — For p = 1, the above equation reads: 



92 



d 



_+4[2(n-K)-?2]^+4,- 



logn„(<r; k) 



+6 ( ^logn„(?;K) 



8nK. 



This can be recognised as the Chazy I equation (see Appendix E) 

if'" + 6{ip'f + 4 [2(n -k)- ip' + Ac;ip - Snn = 0, 



where 



d 

<P('?) = ^logn„ (<?;k). 



(4.29) 

(4.30) 
(4.31) 



Equation (4.30) can further be reduced to the fourth Painleve equation in the Jimbo- 
Miwa-Okamoto a form (Forrester and Witte 2001, Tracy and Widom 1994): 

Piv: {if")'^ -A{ip-^(p')^ +'i(p'{<p' + 2n){<f' -2k) = 0, (4.32) 

see Appendix E for more details. The boimdary condition to be imposed at infinity is 

^(0 ^!^(l + 0(,-i)). (4.33) 

!-oo <r 

Equations (4.20), (4.21), (4.28) and their one-point reductions Eqs. (4.22), (4.23), (4.31) 
and (4.32) are the main results of this subsection. They will play a central role in the 
forthcoming analysis of the replica approach to GUE. 

I Indeed, the boundary conditions at infinity, 

n„|p(?;re) 

do distinguish between the correlation functions characterised by different k's, see Eq. (4.1). 
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4-2. Laguerre Unitary Ensemble (LUE) 

The correlation function of characteristic polynomials in LUE is defined by the formula 
nn|p('i;«) = i / n MA,e-A.AJ n(?„-A,)- ) •A2(A), (4.34) 
where 

n n 

= / n (^^j- e"'^^') • = n + mj + ^) (4.35) 

is the normalisation constant, and it is assumed that v > —1. The associated r function 
equals 

r„^(^,K;t) = l / n(dA,e-^^+''(*=^^)A,^n(?„-A,)'^<^yA^(A). (4.36) 

% j = l V a=l / 

In the above definitions, the superscript L stands for LUE but it will be omitted from 
now on. 

4-2.1. Virasoro constraints 

In the notation of Section 3, the definition Eq. (4.34) implies that || 

/(A) = l ^ ak=6k,i, (4.37) 

g{X) = X-i/ I-)- bk = -iySk,o + 5k,i, (4.38) 

D = R+ h-s. dim(c') = 0. (4.39) 

This brings the following Virasoro constraints Eqs. (2.15) - (2.21) for the r function 
Eq. (4.36): 



dtq dtq^i ^ dtq-m 



T„{<;, k; t) = 'hg T„{<;, k; t). (4.40) 



where 23, is defined by Eq. (4.8) and jtq{t) is the Virasoro operator given by Eq. (2.21). 

In what follows, we need the three lowest Virasoro constraints for g = 0, = 1 
and q = 2. Written for logT„(<j, k; t), they read: 

^jtj-^ - -ET \ logr„(^,K;;t) +n(n + i/ + At) = So log r„(^, «;;*), (4.41) 



K3' 



'dtj dti 



+nt?i(^,K) =Si logT„(^,K;t), (4.42) 



dtx ^ dtl 

2 



d 

logr„(^,K;;t) ) + n t?2(^, «) = logr„(^, k; t). (4.43) 



,dti 

Notice that dim(c') = follows from Eq. (3.70) in which Zo = {0} 
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4-2.2. Toda Lattice hierarchy 



To generate the Toda Lattice hierarchy for the correlation function U^^^{t;;K) 
[Eq. (4.34)] of characteristic polynomials we apply the projection formula 

nn|p(^;«) = ^^n('i,'«;<)[^^, (4-44) 



in which the r function is defined by Eq. (4.36), to the first and second equation of the 
t-dependent TL hierarchy: 

TT 1.. m T„+i(t)r„-i(t) 

't„+i(<)' 



TLo : 



d , ..X T„+i(t)T„_i *) d 

-— r— logr„(t) = 277^ ^ 

dtidt2 T^{t) dti 



Tn-l{t) 



(4.45) 
(4.46) 



see Eqs. (4.15) and (4.16). 



(i) To derive the first equation of the TL hierarchy for n^|p(<j; k) from Eqs. (4.44) 
and (4.45), we have to determine 



92 
dtj 



with the help of the Virasoro constraints. Differentiating Eq. (4.41) over ti and setting 
t = afterwards, we obtain: 

52 



dtf 



logr„(t) = (1 - So) J- log r„(t) 



t=o 



dti 



t=o 



Second, we set t = in Eq. (4.41) to identify the relation 

d 



dti 



log Tn{t) 



t=0 



n {n + u + k) - "Bo logT„(0). 



Combining Eqs. (4.47) and (4.48), we conclude that 

52 
dti 



^logr„(*) 



t=o 



nin + u + K)-'Boil-'Bo) logT„(0). 



(4.47) 



(4.48) 



(4.49) 



Finally, substituting this result back to Eq. (4.45), and taking into account Eqs. (4.44) 
and (4.35), we end up with the first TL equation 



L 

TLi : 



®o(®o-l) logn„|j,(^;K) 



n{n + u) 



n2|^(^;K) 



1 



riK 



(4.50) 



written in the space of physical parameters 

Notice that the above equation becomes more symmetric if written for the correlation 
function 

n„ip(<?) = n„|p(^) n (4.51) 



a=l 



The corresponding TL equation reads: 

~L 

TLi : So(®o-l) log n„|p(.i; k) = n(n + !/) 



nn+i|p(^;K) n„_i|p(«^;K) 

n2|^(^;K) 



1 . (4.52) 
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(ii) The second equation of the TL hierarchy for n„|p(<;; k) can be derived along the 

same hncs. Equation (4.46) suggests that, in addition to the derivative d/dtilogTn at 
t = given by Eq. (4.48), one needs to know the mixed derivative 



dtidt. 



■logr„(t) 



It can be calculated by combining Eq. (4.41) differentiated over t2 with Eqs. (4.42) and 
(4.48). Straightforward calculations bring 



■ log T„(t) 



. =(2-^o)^logr„(t) 

t=0 Ot2 



t=0 



where 

^logr„(t) 



t=o 



dtidt2 

= n [(2n + v + K){n + v + K)+ i?i(<?, k)] 

{2n + v + k) 3o + ^il logT„(0). 



(4.53) 



(4.54) 



The final result reads: 

(So - 2) [(2n + 1/ + k) So + ■Bil logn„|p(^;K) 



TL, 



= n{n + u) 



n„+i|j,(^; /«) n„_i|p(^; k) 



2 (2n + 1/ + k) - So log 



n„+i|p(^;K) 



.n„_i|p(^;K)^ 

-2n (2n + 1/ + k) (n + z/ + k) -n?9i(^, k). (4.55) 

This equation takes a more compact form if written for the correlation function n„|p 
defined by Eq. (4.51): 

~L 

TL2 : (So - 2) [(2n + + k) So + SiJ log n„|p(^; k) 
n„+i|p(?;/«) n„_i|p(^;/c) 



n(n + v) ■ 



n2 {<;-k) 



2(2n + !/) - So log 



nn+i|p(^; /«) 



-2n{n + u) (2n + !/ + «). 



(4.56) 



4-2.3. KP hierarchy and Painleve V equation 



The same technology is at work for projecting the KP hirerachy Eq. (2.8) onto t = 0. 
Below, only the first KP equation 

2 

= (4.57) 



52 ^ 92 

4- 



log Tn{t) + 6 [ ^ log T„(t) 



will be treated. Notice that no superscript (s) appears in Eq. (4.57) as the LUE 
confinement potential does not depend on the matrix size n. To make the forthcoming 
calculation more efficient, it is beneficial to introduce the notation 




logT„(i) 



, T = logT„(0), 



(4.58) 



t=o 
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which brings the KP equation Eq. (4.57) projected onto t = to the form 
Tun + 3T22 - 4Ti3 + ST^ = 0. 



(4.59) 



(i) First, wc observe that Tn and Tun can be determined from the following chain of 
relations, obtained by repeated differentiation of the first Virasoro constraint Eq. (4.41) 
with respect to ti: 



Hence, 



and 



Ti = n (n + z/ + k) - So T, 
Tn =(1-:Bo)Ti, 
Tin = (2 — !Bo) Tn, 
Tun = (3 — !Bo) Tm. 



Tn=n{n + v + K)-{l- 'Bo)'Bo T 

Tnn = 3\n{n + u + K)-{3- Bo)(2 - ^o)(l - ^o)^o T. 



(4.60) 

(4.61) 
(4.62) 



(ii) Second, to determine T13, we differentiate the first Virasoro constraint Eq. (4.41) 
with respect to t^, and make use of the second and third constraints as they stand 
[Eqs. (4.42) and (4.43)] to obtain: 

ri3 = (3-So)r3, 

T3 =(2n + i/ + K)r2+ i9i(^,K)Ti+T2+Tn-®2T + n^?2(<i,K), (4.63) 

T2 = (2n + i/ + K)Ti -!BiT + ni9i(^,K). 

Although easy to derive, an explicit expression for T13 is too cumbersome to be explic- 
itly stated here. 

(iii) To calculate T22, the last unknown ingredient of Eq. (4.59), we differentiate the 
first and second Virasoro constraints [Eqs. (4.41) and (4.42)] with respect to t2 to 
realize that 



J r22 = 2T3 + {2n + u + K) T12 - Bi T2, 

\ri2 = (2-So)T2, 

Combining Eqs. (4.63) and (4.64), one readily derives a closed expression for r22 



(4.64) 



Finally, wc substitute so determined Tnn, Tn, ri3 and r22 into Eq. (4.59) to generate 
a nonlinear differential equation for logn„|p(<j; k) in the form 



KP 



1 • 



+(2n + iy + k)3i(2So - 1) - 2:B2(2!Bo + 1) 



logn„|p(^;«;) 



+6 logn„|p(.i; k))' - 4 (So logn„|p(<;; a^)) logn„|p(^; k)) 

= n [(1/ + k) i?i (^, k) + t?2 (<? , k)] . (4.65) 
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Remark. — For p = 1, the above equation reads: 



+ {2n + u- k) 2^ 



9^2 



d_ 



d 



logn„|p(c; k) 

g^ + ^^d^ ) logn„|p(?;K) 



It can further be simplified if written for the function 

d 

= logn„(c; k) - riK. 



(4.66) 
(4.67) 



^21 ^' 



Straightforward calculations yield: 
ci'^tp'" + Q,^" - [c^ - 2(2n + - k) c + 4«;n + (z/ + k)^ 

+ [2(2«; - n) - 1/ - 2^] + 6 ^(1^')^ - ^w' = 2n«;(n + i/). (4.68) 
This can be recognized as the Chazy I form (see Appendix E) of the fifth Painleve 
transcendent. Equivalently, satisfies the Painleve V equation in the Jimbo-Miwa- 
Okamoto form (Forrester and Witte 2002, Tracy and Widom 1994): 

Pv : {S'p"f - - <i^' + '^W? + (2n + - k)^'] ' 

+4<^'((^' + n)^ + n + z/)(v>' - k) = 0. 
Both equations have to be supplemented by the boundary condition 



(4.69) 



(4.70) 



following from Eq. (4.67) and the asymptotic analysis of Eq. (4.34). Equations (4.50), 
(4.52), (4.55), (4.56), (4.65) and (4.69) represent the main results of this subsection. 

J^.S. Discussion 

This Section concludes the detailed exposition of integrable theory of correlation 
functions of RMT characteristic polynomials. Among the main results derived are: 

• The multivariate first [Eqs. (4.20) and (4.52)] and second [Eqs. (4.21) and (4.56)] 
equations of the Toda Lattice hierarchy ^ which establish nonlinear differential 
recurrence relations between "nearest neighbor" correlation functions n„|p(^;K) 

and n„±i|p(i;; /«), and 

• The nonlinear multivariate differential equations [Eqs. (4.28) and (4.65)] satisfied 
by n„|p(^; k) alone. These can be considered as multivariate generalisations of the 
corresponding Painleve equations arising in the one-point setup p = 1 [Eqs. (4.32) 

and (4.69)] 

Other nonlinear multivariate relations between the correlation functions n„|p(^; k) and 
n„-|-q|p(ij; k) can readily be obtained from the modified Toda Lattice and Kadomtsev- 
Petviashvili hierarchies listed in Section 3.4. 

Finally, let us stress that a similar calculational framework applies to other /3 = 2 
matrix integrals depending on one (Osipov and Kanzieper 2007) or more (Osipov and 
Kanzieper 2009; Osipov, Sommcrs and Zyczkowski 2010) parameters. The reader is 
referred to the above papers for further details. 



^ See also their single Vciriable reductions Eqs. (4.22) and (4.23) derived for the GUE. 
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5. Integrability of Zero-Dimensional Replica Field Theories 

5.1. Introduction 

In this Section, the integrable theory of CFCP will be utilised to present a tutorial 
exposition of the exact approach to zero dimensional replica field theories formulated 
in a scries of publications (Kanzicpcr 2002, Splittorff and Verbaarschot 2003, Osipov 
and Kanzieper 2007). Focussing, for definiteness, on the calculation of the finite- 
average eigenlevel density in the GUE (whose exact form (Mehta 2004) 

^^^'^ = 2A^r(!v)V^ '^'' [H'^mN-iie) - HNie)H'^_,ie)] 

has been known for decades), we shall put a special emphasis on a comparative analysis 
of three alternative formulations - fermionic, bosonic and supersymmetric - of the 
replica method. This will allow us to meticulously analyse the fermionic-bosonic 
factorisation phenomenon of RMT spectral correlation functions in the fermionic and 
bosonic variations of the replica method, where its existence is not self-evident, to say 
the least. 

5.2. Density of eigenlevels in finite-N GUE 

To determine the mean density of eigenlevels in the GUE, we define the average one- 
point Green function 

G(.;7V) = (tr(.-JC)-^)„^^^^^ (5.1) 
that can be restored from the replica partition function (n e K"*") 

ZW(.; N) = (det±"(. - ?C))„^^^^^ (5.2) 
through the replica limit 

Equation (5.2) can routinely be mapped onto cither fermionic or bosonic replica field 
theories, the result being (see, e.g., Kanzieper 2010) 



{z; N) = - Z- / (D„Q) e-*'-''^' det^ {iz - Q) (5.4) 

Cn J 



and 



4-)(^;7V) = l /(2)„0)e-*--''°'det-^(^-0). (5.5) 

Cn J 

Both integrals run over nxn Hermitean matrix O; the normalisation constant c„ equals 

Cn = Ji^nQ) e-''-^\ (5.6) 

By derivation, the replica parameter n in Eqs. (5.4) and (5.5) is restricted to integers, 
n e Z+. Notably, Eqs. (5.4) and (5.5) are particular cases of the correlation function 
n„|p(^; k) studied in previous sections. 

5.2.1. Fermionic replicas 

Indeed, comparison of Eq. (5.4) with the definition Eq. (4.1) yields 

Zi+)(z;iV)- (-Z)"^n„(iz;7V), (5.7) 

where the shorthand notation n„(2;;iV) is used to denote n^|^(2;; A''), in accordance 
with the earlier notation in Eqs. (4.22) and (4.23). This observation results in the 
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Painleve IV representation of the fermionic replica partition function [see Eqs. (4.31) 
and (4.32)]: 

^logZi+\z;N) = i^{t;n,N) , (5.8) 

where ip{t\ n, N) is the fourth Painleve transcendent satisfying the equation 

{ip"f - 4{ip -tip'f + 4ip'{ip' + 2n)((/?' -2N) = (5.9) 

subject to the boundary conditions + 

nN 

ip{t;n,N) —, \t\^oo, tGC. (5.10) 

Here and above, n G Z_|_. 

Equations (5.8) and (5.9) open the way for calculating the average Green function 
G{z; N) via the fermionic replica limit 

G(z]N)=lim--^Zi+Hz]N) = Llim-ip(t]n,N) . (5.11) 

For the prescription Eq. (5.11) to be operational, the Painleve representation of 
7}h^{z\N) should hold * for n G Notice that for generic real n, the fermionic 

replica partition function zi'^^(2;; N) is no longer an analytic function of z and exhibits 
a discontinuity across the real axis. For this reason, the Painleve equation Eq. (5.9) 
should be solved separately for £flet < {^mz > 0) and Diet > (3m 2; < 0). 

Replica limit and the Hamiltonian formalism. — To implement the replica limit, we 
employ the Hamiltonian formulation of the Painleve IV (Noumi 2004, Forrester and 
Witte 2001) which associates ip{t; n, N) with the polynomial Hamiltonian (Okamoto 
1980a) 

^{t; n, N) = Hi {P, Q, t} = {2P + Q + 2t)PQ + 2nP - NQ (5.12) 

of a dynamical system {Q,P,Hf}, where Q = Q{t;n,N) and P — P{t;n,N) are 

canonical coordinate and momentum. For such a system, Hamilton's equations of 
motion read: 

Q= + ^ = Q(Q + 4P + 2i) + 2n, (5.13) 
oP 

P= - ^ = + 2P + 2t) + N. (5.14) 

Since 

G{z;N) = LVLm.-Hi{P,Q,t} , (5.15) 

n^O n t=iz 

we need to develop a small-n expansion for the Hamiltonian Hf {P,Q,t}. Restricting 
ourselves to the linear in n terms, 

Hi {P, Q, t} = nH^it; N) + 0{n^) (5.16) 

and 

P{t:n, N) = pa{t; N) + npi{t; N) + 0{n^), (5.17) 

Q{t; n, N) = qo{t; N) + nqi{t; N) + 0{n^), (5.18) 

+ Equation (5.10) follows from Eqs. (5.8), (5.7) and the footnote below Eq. (4.28). 

* Previous studies (Kanzieper 2002, Osipov and Kanzieper 2007) suggest that this is indeed the case. 
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we conclude that q(){t:N) = 0. This derives directly from the expansion Eq. (5.16) in 
which absence of the term of order 0{nP) is guaranteed by the normalisation condition 

Zi'^\z; N) = l. As the result f, 

G{z;N) = LH[^\Lz;N), 

where 



(5.19) 



Hx{t- N) = 2pogi(po +t) + 2po - Nqi, 
Hiit;N)=2poqi. 

Here, po = Po{t; N) and qi = qi{t; N) are solutions to the system of coupled first 
equations: 

Pa = -2pl- 2pat + N, 
qi = 4po?i + 2qit + 2. 

Since the initial conditions are known for Hi {t; N) rather than for po{t; N) and qi 
separately, below we determine these two hmctions up to integration constants. 
The function po{t; N) satisfies the Riccati differential equation whose solution 



t 



5.20) 
5.21) 
order 

5.22) 

t;N) 
is 

5.23) 

5.24) 

5.25) 
5.24) 

will simplify the formulae to follow. Second, the solution Eq. (5.23) for po{t; N) actually 
depends on a single constant (either Ci/c2 or C2/C1) as it must be. 

To determine qi{t; N), we substitute Eq. (5.23) into the second formula of Eq. 
to derive: 

dt 



u+{t) 



where 



u+{t) = ci£>_jv_i(i\/2) + C2i-L)^DN{itV2) 
is, in turn, a solution to the equation of parabolic cylinder 

il+it) -{2N + 1+ t^)u+it) = 0. 
Two remarks are in order. First, factoring out (— i)^ in the second term in Eq. 



qi{t;N)=2uiit) 



Making use of the integration formula (sec Appendix F) 

dt _ 1 aiD_N-iitV2) + a2{- L)'^DN{itV2) 



u\{t) V2 u+{t) 
where two constants ai and a2 are subject to the constraint 

Cia2 — C2ai = 1, 
we further reduce Eq. (5.26) to 

qiit;N) = V2u+{t) \aiD_N-i{tV2) + a2{-L)^ DN(itV2) 

Equations (5.21), (5.23), (5.29) and the identity 

u+{t) -tu+{t) = -V2 \ciD_N{ty/2) - C2{-l)^~'^NDN-\{it\/2) 



5.22) 
5.26) 



5.27) 

5.28) 
5.29) 

5.30) 



(obtained from Eq. (5.24) with the help of Eqs. (F.7) and (F.8)) yield Hi{t;N) in the 
form 



Hi{t;N) = -2 [ai£'_jv-i(t\/2) + Q!2(-i)^^Jv(^i^/2)] 

X [ciD_NitV2) - C2i-if-^NDN-i{LtV2)] . 
t We will drop the superscript (f) wherever this does not cause a notational confusion. 



(5.31) 
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Notice that appearance of four integration constants (ci, C2, ai and in Eq. (5.31) is 
somewhat illusive: a little thought shows that there is a pair of independent constants, 
either (ci/c2, Q:2C2) or their derivatives. 

To determine the imknown constants in Eq. (5.31), wc make use of the asymptotic 
formulae for the functions of parabolic cylinder (collected in Appendix F) in an at- 
tempt to meet the boundary conditions | 

H,{t;N)^j, H,{t;N)r^-^, 1*1^00, t G C. (5.32) 

Following the discussion next to Eq. (5.11), the two cases Diet < and Diet > will 
be treated separately. 

• The case ^Ret <0. Asymptotic analysis of Eq. (5.31) at t — )• — oo yields 



57 = (-" ivr- 

so that 

Hi{t; N) = 2D_N-i{-tV2) [aiCi£>_jv(-t\/2) - i'^ N D n -i{LtV2) \ . (5.33) 

Here, wc have used Eq. (F.5). To determine the remaining constant aici, wc 
make use of the boundary condition Eq. (5.32) for t — >■ ±ioo — 0. Straightforward 
calculations bring aiC\ = 0. We then conclude that 

Hi{t;N) = -2{-i)^-'^ND_N-i{-t'^)DN-ii-itV2), Dlet<0. (5.34) 



The case ^Ret > 0. Asymptotic analysis of Eq. (5.31) at t ^ +00 yields a2 = so 
that 



Hi{t; N) = 2D_N-i{tV2) 



Cl 



D_N(tV2) - (-i)^-'^NDN-iUt\/2) 

C2 



(5.35) 



To determine the remaining constant Ci/c2, we make use of the boundary condition 
Eq. (5.32) for t ±ioo + 0. Straightforward calculations bring Ci/c2 = 0. We 
then conclude that 

Hi{t;N) = -2{-L)^-^ND_N-i{tV2)DN-i{itV2), Diet > 0. (5.36) 

The calculation of Hi{t: N) can be summarised in a single formula 

Hi{t; N) = -2{-if-^ND_N-i{<JittV2)DN-i{LautV2), (5.37) 

where = sgn (it) = sgn Die t denotes the sign of Die t. In terms of canonical 
variables po{t;N) and qi{t;N), this result translates to 

qi{t:N)= -V2au{-LfD-NMtaitV2)DN{itauV2). (5.39) 

Now Hi{t;N) can readily be restored by integrating Eq. (5.37). We proceed in three 
steps, (i) First, we make use of differential recurrence relations Eqs. (F.7) and (F.8) 
and the Wronskian Eq. (F.4) to prove the identity 

iND-N-i{(Titt\/2)DN-iiiaittV2) = - D-N{(nttV2)DN{L(TittV2). (5.40) 
I Equation (5.32) is straightforward to derive from Eqs. (5.16), (5.12) and (5.10). 
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The latter allows us to write down Hi{t; N) as 

Hi{t; N) = -2 + 2{-i)^ D_N{cTittV2)DN{i<JutV2). (5.41) 
(ii) Second, it is beneficial to employ the differential equation Eq. (F.3) to derive 



■^D-N{(TittV2) = {2N-1+ e)D-N{(JittV2), 



DN{i(Titty/2) ={2N+1 + t'^)DN{iaittV2). 



(5.42) 
(5.43) 



These two relations imply 

2D_N{aittV2)DN{LautV2) = D_N{(TufV2) — DN[LautV2 



-Dn {iaittV2) -^D-N {(TutV2) 



(5.44) 



so that 

Hi{t;N) = -2 



D_N{(TittV2)-^DN{iautV2) - DN{L(TittV2)-^D_N{(TittV2) 



(iii) Third, we integrate the above equation to obtain 

Hi{t;N) = -2t+{-i)^Wt \D_N{(TutV2),DN{iaittV2) 



(5.45) 



(5.46) 



Here, the integration constant was set to zero in order to meet the boundary conditions 
Eq. (5.32) at infinities. The notation stands for the Wronskian 



(5.47) 



Average Green function and eigenlevel density. — Now, the average one-point Green 
function readily follows from Eq. (5.19): 



G{z; N) = 2z + (-O^W, L»_Ar(-izcr^v^), DN{za^V2) 
Here, = sgnJmz denotes the sign of imaginary part of z = —it. 



(5.48) 



The average density of eigenlevels can be restored from Eq. (5.48) and the relation 



Q^{e) = -3mG{e + la^Q;N). 

TT 

Indeed, noticing from Eqs. (F.2) and (5.48) that 

am [{-irD.,,{-iea.V2)\ = "^i^^^^ 
we conclude, with the help of Eq. (F.l), that 



/2 / dTT- -e 



-N-l. 



(5.49) 



(5.50) 



am 



(-i)'^-D-jv(-ieCT^\/2)] 



,N-l 



'/^HN-i{e). (5.51) 



2^/2r(7V) 

Here, HM-i{e) is the Hermite polynomial appearing by virtue of the relation 
DN{zy/2) = e -^NpT' 



(5.52) 
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Consequently, 



(5.53) 



2^r(7V) 

Taken together with Eqs. (5.49) and (5.48), this equation yields the finite-A'' average 
density of eigenlevels in the GUE: 
1 



£'Jv(e) 



2^r(Ar)^ 



1 



(5.54) 



While this result, obtained via the fermionic replica limit, is seen to coincide with the 
celebrated finite- AT formula (Mehta 2004) 

1 



[H'^{e)HN^,{e) - ifjv(e)i?^_i(6)] (5.55) 



2^r(Ar)v^ 

originally derived within the orthogonal polynomial technique, the factorisation 
phenomenon (as defined in Section 1.2) has not been immediately detected throughout 
the calculation of either G{z; N) or qn{c)- We shall return to this point in Section 5.3. 

5.2.2. Bosonic replicas 

Comparing Eq. (5.5) with the definition Eq. (4.1), we conclude that 

Zi-\z;N)=Iln{z;-N), (5.56) 

where 3mz ^ 0. The shorthand notation ]In{z;—N) is used to denote n^|^(2;; — A''), 
in accordance with the earlier notation in Eqs. (4.22) and (4.23). Consequently, 
the Painleve IV representation of the bosonic replica partition function reads [see 
Eqs. (4.31) and (4.32)]: 

^;log4-)(^;iV) = <^(i;n,-Ar) 



dz 



t=z 



(5.57) 



where ijj{t;n,N) = ip{t;n,—N) is the fourth Painleve transcendent satisfying the 
equation 

~ - ttp'f + Aijj'i^' + 2n)(V'' + 2N) = 
subject to the boundary conditions 

nN 



^{t;n,N) 



t 



\t\ DO, 



ieC\] 



(5.58) 
(5.59) 



Here and above, n e 



The average Green function G{z; N) we are aimed at is given by the bosonic replica 
limit 



G{z; N) = - lim ^ ^4"^ (^i N) = - Urn ^ ^{t; n, N) 



n->-0 n 



t=z 



(5.60) 



To implement it, we assume that the Painleve representation of Zn \z;N) holds for 
n e M+. 

Replica limit and the Hamiltonian formalism. — Similarly to our treatment of the 
fermionic case, we employ the Hamiltonian formulation of the Painleve IV (Noumi 2004, 
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Forrester and Witte 2001) which associates tp{t; n, N) with the polynomial Hamiltonian 
(Okamoto 1980a) 

^P{t; n, N) = H^{P, Q, t} = {2P + Q + 2t)PQ + 2nP + NQ (5.61) 

of a dynamical system {(5,P, TJb}, where Q = Q(t-,n,N) and P — P{t;n,N) are 
canonical coordinate and momentum. For such a system, Hamilton's equations of 
motion read: 

Q = + -gp- = Q{Q + + '2t) + 2n, (5.62) 

P = - ^ = -P{2Q + 2P + 2t) - N. (5.63) 

Owing to Eq. (5.60), we need to develop a small-n expansion for the Hamiltonian 
i?b{Q,P,i}: 

H^{P,Q,t} = nH^^\t-N) + 0(n2). (5.64) 
Being consistent with yet another expansion 

P{t;n,N) = pQ{t;N) + npx{t;N) + 0{n^), (5.65) 

Q{t;n,N)=nqi{t;N) + 0{n'^), (5.66) 
it results in the relation f 

G{z;N) = -H^^\z;N), (5.67) 

where 

Hi{t- N) = 2poqiipo +t) + 2po + Nqi, (5.68) 

Hi{t;N)=2poqi. (5.69) 

Here, po = Po{t; N) and qi — qi{t; N) are solutions to the system of coupled first order 
equations: 

qi = Apoqi+2qit + 2. 

Since the initial conditions are known for Hi{t; N), rather than for po{t; N) and qi{t; N) 
separately, below we determine these two functions up to integration constants. 
The function po{t; N) satisfies the Riccati differential equation whose solution is 



U-{t) 



(5.71) 



where 

u-{t) = cii^DN-i{tV2) + C2i^D_N{it\f2) (5.72) 

is, in turn, a solution to the equation of parabolic cylinder 

u_ {t) + (2iV - 1 - t^)u- {t) = 0. (5.73) 

Factoring out in the second term in Eq. (5.72) will simplify the formulae to follow. 

To determine qi{t; N), we substitute Eq. (5.71) into the second formula of Eq. (5.70) 
to derive: 

q,{t;N)=2ulit)J Jj-^. (5.74) 

Making use of the integration formula (see Appendix F) 

dt _ I aiDN-iitV2) +a2D_N{itV2) 



I 



t We will drop the superscript (b) wherever this does not cause a notational confusion. 



(5.75) 



5 Integrability of Zero-Dimensional Replica Field Theories 



39 



where two constants ai and a2 are subject to the constraint 

Cia2 — C2ai = 1, 
we further reduce Eq. (5.74) to 

qi{t;N) = V2u-{t) aiDN-i{t^) + a2D_N{itV2) 

Equations (5.69), (5.71), (5.77) and the identity 

u-{t) -tu-{t) = -y/2 \cil^DN{tV2) + C2L^+^ND_N-i{ltV2) 



(5.76) 
(5.77) 

(5.78) 



(obtained from Eq. (5.72) with the help of Eqs. (F.7) and (F.8)) yield Hi{t;N) in the 
form 



(5.79) 



To determine the unknown constants in Eq. (5.79), we make use of the asymptotic 
formulae for the functions of parabolic cylinder (collected in Appendix F) to satisfy 
the boundary conditions [see Eq. (5.59)] 



N ■ N 



|i|^-oo, ieC\] 
The two cases drat < and 3m t > should be treated separately. 



(5.80) 



The case 3m.t < 0. Asymptotic analysis of Eq. (5.79) at t ^ —too yields ci = 0, 
so that 



Hi{t;N) = 2L^+^ND_N-i{itV2) [£)jv-i(t\/2) - a2C2D_N {its/2) 



(5.81) 



To determine the remaining constant a2C2, we make use of the boundary condition 
Eq. (5.80) for t — >■ ±oo — iO. Straightforward calculations bring a2C2 = 0. We then 
conclude that 



Hi{t]N) = 2 L^+'^N D_N-i{itV2)DN-i{tV2), Omt < 0. 



(5.82) 



The case 3mt > 0. Asymptotic analysis of Eq. (5.79) at t — > +ioo yields 



— = - 

C2 



:\N-1 



27r 



(AT-l)! 



(5.83) 



so that 



Hi{t; N) = -N\J-D_N-ii-itV2) 



D_N{itV2) + aiCi{-L) 



N- 



-i(N-l)l 



27r 



D-N{-LtV2) 



(5.84) 



To determine the remaining constant ctiCi, we make use of the boundary condition 
Eq. (5.80) for t — >■ ±oo + iO. Straightforward calculations bring 



aici = (-i) 



;\N-1 



2-K 



(5.85) 



{N-l)\ 
We then conclude that 

Hi{t;N) = 2i^+'^N D-N-i{-LtV2)DN-i{-tV2), 3mt>Q. (5.86) 
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The calculation of Hi {t; N) can be summarised in a single formula 

Hi{t;N) = -2L^-^ND_N-i{-''t(TtV2)DN-i{-t(TtV2), (5.87) 

where crt = sgn!llmi denotes the sign of 3xnt. In terms of canonical variables PQ{t; N) 
and qi{t;N), this result translates to 

qi{t;N) = V2ati^ D_N{-itatV2)DN-i{-tatV2). (5.89) 
In view of Eq. (5.67), the latter result is equivalent to the statement 

— Giz;N) = ~H[^\t;N) ^ 2 L^-^ND_N-i{-iza,V2)DN-ii-za,V2). (5.90) 

OZ t—z 

This expression, obtained within the bosonic replicas, must be compared with its 
counterpart derived via the fermionic replicas [Eqs. (5.19) and (5.37)]: 

— G{z;N) = -H[^\t;N) =2{-L)''-^ND_N-i{-iza,V2)DN-i{za,V2). (5.91) 

OZ t=iz 

As the two expressions coincide, we are led to conclude that the bosonic version of 
the replica limit reproduces correct finite- A/' results for the average Green function and 
the average density of cigcnlcvcls as given by Eqs. (5.48) and (5.55), respectively. 
Again, as is the case of a fermionic calculation carried out in Section 5.2.1, the 
factorisation property did not show up explicitly in the above bosonic calculation. 
We defer discussing this point until Section 5.3. 



5.2.3. Supersymmetric replicas 



The very same integrable theory of characteristic polynomials is at work for a "su- 
persymmetric" variation of replicas invented by Splittorff and Verbaarschot (2003). 
These authors suggested that the fermionic and bosonic replica partition functions 
(satisfying the fermionic and bosonic Toda Lattice equations f, respectively) can be 
seen as two different branches of a single, graded Toda Lattice equation. Below we 
show that the above statement, considered in the context of GUE, is also valid beyond 
the first equation of the Toda Lattice hierarchy. 

First (graded) TL equation. — Equations (5.7) and (4.22) imply that the fermionic 
replica partition function zl^\z;N) satisfies the first TL equation in the form: 

^log2;(-)(.;Ar) = -2n ^^^.^^ 1 

Together with the initial conditions Zq^\z;N) = 1 and 
Z^^\z;N) = {-i)^Ilf,^{iz;N) = ^ [ d\e-^\z + iX)^ 

= 2-'''^e''"/^DM{zV2), (5.93) 

this equation uniquely determines fermionic replica partition functions of any order 
(n > 2). Here, Dj^ is the function of parabolic cylinder of positive order (see 
Appendix F). 




t Notice that SplittorflF and Verbaarschot (2003) use the term "Toda Lattice equation" for the first 
equation of the TL hierarchy. 
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The first TL equation for tlie bosonic replica partition function Zn \z;N) follows 
from Eqs. (5.56) and (4.22), 

^logZi )(.;iV) = +2n(^ ^^^^ ij . (5.94) 

Together with the initial conditions Zq~\z; N) = 1 and 

z[-\z;N) = U'^],iz;-N) = ^ f dXe-^\z-X)-^ 

VT^ Jr 

= {-La,)^2^/^e-''^^ D-Nhiza,V2), (5.95) 

where = sgnJmz denotes the sign of Jmz, this equation uniquely determines 
bosonic replica partition functions of any order (n > 2). Here, D-jv is the function of 

parabolic cylinder of negative order (see Appendix F). 
Further, defining the graded replica partition function as 

(z\-\z;N), neZ- 

Zn{z;N) = i 1, n = (5.96) 

[z(J(z;iV), n€Z+, 

we spot from Eqs. (5.92) and (5.94) that it satisfies the single (graded) TL equation 

^logM-;A^) = -2n(^ ^^j^^ Ij. (5.97) 

Here, the index n is an arbitrary integer, be it positive or negative. The first graded 
TL equation must be supplemented by two initial conditions given by Eqs. (5.93) and 
(5.95). 

Second (graded) TL equation. — Equations (5.7), (5.56) and (4.23) imply that both 
fermionic and bosonic replica partition functions zl^\z;N) additionally satisfy the 
second TL equation 

(l - z^) ^ log Zi^Hz- N) = n ^'nM-^^)^'nMz;N) d_ ( 0£i^\ 
'dz)dz''^^- ^'^""^ zi^^\z;N) Oz''^[z^^)^^,,N))- 

(5.98) 

Consequently, the graded replica partition function Zn{z; N) defined by Eq. (5.96) is 
determined by the second graded TL equation 

1 ^ ^ ^ , 7 / ^^^ ^n+iiz;N) Zn-i{z;N) d f Z^+^jz; N) 

'-'d-z)d-z = " — w^) — ^ [z;:^^), 

(5.99) 

supplemented by two initial conditions Eqs. (5.93) and (5.95). 

Replica limit of graded TL equations. — To determine the one-point Green function 
G{z; N) within the framework of supersymmetric replicas, one has to consider the 
replica limit 

G(z; N) = lim ZJz; N). (5.100) 
n->0 n oz 

The first and second graded TL equations bring 

G'{z; N) = 2-2 Z-i{z; N) Zi{z; N), (5.101) 
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and 

G{z; N) = zG'{z; N) + Z--i{z; N)Z[{z; N) - Z'_i{z; N)Zi{z; N), (5.102) 

respectively. Combining the two equations, we derive 

G{z;N)=2z-2zZ-iZi+W^[Z-i, Zi], (5.103) 

where Wz is the Wronskian defined in Eq. (5.47); the prime ' stands for the derivative 
d/dz. Interestingly, the second graded TL equation has allowed us to integrate 
Eq. (5.102) at once! 

The resulting Eq. (5.103) is remarkable: it shows that the average Green function 
can solely be expressed in terms of bosonic Z_i(z;iV) and fermionic Zi{z;N) replica 
partition functions with only one flavor. This structural phenomenon known as the 
'factorisation property' was first observed by Splittorff and Verbaarschot (2003) in the 
context of the GUE density-density correlation function. Striking at first sight, the 
factorisation property appears to be very natural after recognising that fermionic and 
bosonic replica partition functions are the members of a single graded TL hierarchy. 

To make Eq. (5.103) explicit, we utilise Eqs. (5.93) and (5.95) to observe the identity 

Zi] = 2z Z-i Zi + i-ifWz [D_N{-LzazV2), DN{za,V2)^ . (5.104) 
Consequently, 

G{z; N)=2z + \D-N{-iza^V2), £>iv(^(7^\/2)l . (5.105) 

This expression for the average Green function, derived within the framework of 
supersymmetric replicas, coincides with the one obtained separately by means of 
fermionic and bosonic replicas (see, e.g., Eq. (5.48)). Hence, the result Eq. (5.55) 
for the finite- A'' average density of eigenlevels readily follows. 



5. 3. Factorisation property in fermionic and bosonic replicas 

The factorisation property naturally appearing in the supersymmetric variation of the 
replica method suggests that a generic correlation function should contain both compact 
(fermionic) and non-compact (bosonic) contributions. Below, the fermionic-bosonic fac- 
torisation property will separately be discussed in the context of fermionic and bosonic 
replicas where its presence is by far less obvious even though the enlightened reader 
might have anticipated the factorisation property from Eqs. (5.19) and (5.21) for 
fermionic replicas and from Eqs. (5.67) and (5.69) for bosonic replicas. 



Fermionic replicas. -The Hamiltonian formulation of the fourth Painlevc transcen- 
dent employed in Section 5.2.1 is the key. It follows from Eqs. (5.19) and (5.21) that 
the average Green function G{z; N) is expressed in terms of canonical variables po and 
q\ as 



d 

— G{z; N) = -2po{Lz) qi{iz) 



where 



and 



Poii'Z) = 



z + —\ogDN{zV2) 



Qiilz) = - — {-idzfe' 



(5.106) 



(5.107) 



(5.108) 



To derive the last two equations, we have used Eqs. (5.38) and (5.39) in conjunction 
with Eqs. (F.7) and (F.8). 
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With the help of Eq. (5.93), the canonical "momentum" po can be related to the 
fermionic partition function for one flavor, 

Po{iz) = -~logZ[+\z;N). (5.109) 

This is a compact contribution to the average Green function. A non-compact 
contribution is encoded in the canonical "coordinate" qi which can be related to the 
bosonic partition function via Eq. (5.95): 



q,iiz) = -^z[-'\z;N)-^^z[-\z;N), (5.110) 

so that 



This is yet another factorised representation for G'{z;N) [compare to Eq. (5.101)]. 

Bosonic replicas. — To identify both compact and non-compact contributions to the 

average Green function, wc turn to Eqs. (5.67) and (5.69) to represent the derivative 
of the average Green function G{z; N) in terms of canonical variables po and qi as 

^^G{z;N) = -2po{z)qi{z), (5.112) 



where 

(5.113) 



d 
oz 



and 

qi{z) = -^{-ia,)''e-''/^D_t^{-Lza,V2)-^^ [e'''^^Dj,{zV2)] . (5.114) 

To derive the last two equations, we have used Eqs. (5.88) and (5.89) in conjunction 
with Eqs. (F.7) and (F.8). 

With the help of Eq. (5.95), the canonical "momentum" pq can be related to the 
bosonic partition function for one flavor, 

Poiz)^~logZ[-\z;N). (5.115) 

This is a non- compact contribution to the average Green function. A compact 
contribution comes from the canonical "coordinate" qi which can be related to the 
fermionic partition function via Eq. (5.93): 

Qiiz) = ~ z[-\z;N) ^^z[+\z;N), (5.116) 

so that 

^^Giz;N) = ^^^z[+\z;N)^^z[-\z;N), (5.117) 
agreeing with the earlier result Eq. (5.111). 

Brief summary. — The detailed analysis of fermionic and bosonic replica limits per- 
formed in the context of the GUE averaged one-point Green function G{z; N) has 
convincingly demonstrated that the Hamiltonian formulation of the fourth Painleve 
transcendent provides a natural and, perhaps, most adequate language to identify the 
factorisation phenomenon. In particular, we have managed to show that the deriva- 
tive G'{z;N) of the one-point Green function factorises into a product of canonical 
"momentum" 

Po{t; N) = lim P{t; n, N) (5.118) 

n— >-0 
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and canonical "coordinate" 

qJt; N) = lim -Q{t; n, N). (5.119) 

n-5-O n 

As suggested by Eqs. (5.109) and (5.115), the momentum contribution to the average 
Green function is a regular one; it corresponds to a compact contribution in fermionic 
replicas and to a non-compact contribution in bosonic replicas: 

d jz[+\z;N) (fermionic) 
Po ~ 35-log< , . (5.120) 
[Z\~'{z;N) (bosonic) 

On the contrary, the coordinate contribution qi is of a complementary nature: defined 
by a replica-like limit [Eq. (5.119)] it brings in a noncompact contribution in fermionic 
replicas [Eq. (5.110)] and a compact contribution in bosonic replicas [Eq. (5.116)]: 

ff A d iz[-\z;N) (fermionic) 
gi ~ exp / po c^-z X — ^ (5.121) 
\J J 9z [z^+\z;N) (bosonic) 

We close this section by noting that the very same calculational framework should 
be equally effective in performing the replica limit for other random-matrix ensembles 
and/or spectral observables. 



6. Conclusions 



In this paper, we have used the ideas of integrability to formulate a theory of the 
correlation function 

r P 
n„|p(«?; k)= duni'K) n det;:°(fe - ?t) 

•' Q = l 

of characteristic polynomials for invariant non-Gaussian ensembles of Hermitean 
random matrices characterised by the probability measure d/x„(?t) which may well 
depend on the matrix dimensionality n. Contrary to other approaches based on various 
duality relations, onv theory docs not assume "integerness" of replica parameters 
K ~ (Ki,---,Kp) which arc allowed to span k &W. One of the consequences of 
lifting the restriction k G V is that we were unable to represent the CFCP explicitly 
in a closed determinant form; instead, we have shown that the correlation function 
n„|p(<;;/«) satisfies an infinite set of nonlinear differential hierarchically structured 
relations. While such a description is, to a large extent, im,plicit, it does provide 
a useful nonperturbative characterisation of n„|p(^;/t) which turns out to be much 
beneficial for an in-depth analysis of the mathematical foundations of zero-dimensional 
replica field theories. 

With certainty, the replicas is not the only application of a nonperturbative approach 
to CFCP developed in this paper. With minor modifications, its potential readily 
extends to various problems of charge transfer through quantum chaotic structures 
(Osipov and Kanzieper 2008, Osipov and Kanzieper 2009), stochastic theory of density 
matrices (Osipov, Sommers and Zyczkowski 2010), random matrix theory approach to 
QCD physics (Verbaarschot 2010), to name a few. An extension of the above formalism 
to the CFCP of unitary matrices may bring a useful calculational tool for generating 
conjectures on behaviour of the Riemann zeta function at the critical line (Keating and 
Snaith 2000a, 2000b). 

Finally, we wish to mention that an integrable theory of CFCP for ^ = 1 and /3 = 4 
Dyson symmetry classes is very much called for. Building on the insightful work by 
Adler and van Moerbeke (2001), a solution of this challenging problem seems to be 
within the reach. 
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Appendices 

A. Bilineetr Identity in Hirota Form: An Alternative Derivation 

In Section 3.3, the bilinear identity in Hirota form [Eq. (3.27) or, equivalently, 
Eq. (3.34)] was derived from the one in the integral form [Eq. (3.18)]. As will be 
shown below, the latter is not actually necessary. 

An alternative proof of the bilinear identity in Hirota form [Eq. (3.34)] starts with 

the very same Eq. (3.19), 

/ dAr„(A)e''(*'^)e("-i)''(*-*''^)pf)(t;A)Pi")(t';A) 

rfAr„(A) e''(*''^'e"''(*-*':^)pj"^(t; A)Pi,")(t'; A). (A.l) 



IV 

However, instead of employing Cauchy representations [see Eqs. (3.21) and (3.24) in 
the original proof], we assume existence of Taylor expansions 

oo 

^{a-l)v{t-t';X)pM^ti. A) = ^ A^A" (A.2) 

CT=0 

and 

oo 

^av{t-t';X)pin)^^, A) = ^ P<,A^ (A.3) 



CT=0 



to rewrite Eq. (A.l) in the form 



E (^^ Pt\t; A))^ = E 5^ Pk-Kt'; A)}^ . (A.4) 

o-=0 " cr=0 

Here, we have used the scalar product notation defined in Eq. (3.6). The bilinear iden- 
tity is obtained from Eq. (A.4) after expressing all its ingredients in terms of the r 
function Eq. (3.1). 

(i) The coefficients A^- and B„. To determine A„, we employ two identities, 

oo 

g(a-i) vit-t';X) ^J2sk{{a-l){t- t')) (A.5) 

k=0 

and 

P^\t'' A) = ( E ^''rn-k {-[df]) 1 T(r-")(t'), (A.6) 



(rn- 

where [df] denotes 

If) ]-fA. lA 

^ \dt['2dt'^'"'' kdt'^'' 

Equation (A.5) is a consequence of Eq. (2.6) and the definition of Schur polynomials 
(see Table 1). Equation (A.6) follows from the Heine formula 

P^\t';X) = / n (rfA, (A- A,)r„(A,)e''(*'^^^)) • A^(A) (A7) 

and the identity 
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where p^W is an infinite dimensional vector 



PmW = (tr„A, itr„A^ • • • , ^tr„A^ • • • j . (A.9) 
Indeed, substituting Eq. (A. 8) into Eq. (A. 7), we derive: 

, (^-n),,. E / n K- r„(A,) e^'f*'^^^)) . s, (-p„(A)) A^(A) 

m ^ 



fe=0 

A" 



In the last step, we have used the obvious formula 

m m 

(-[St']) n e''^* ''^^ = (-Pm(A)) n ^''^^ (A-ii) 

Having established Eq. (A. 6), we substitute it and Eq. (A. 5) into Eq. (A. 2) to obtain 

A. = E Sfe((a-l)(t-0)«m-.+fc(-[at'])r(r-")(t'). (A.12) 

■^TO (* ) fe=max(0,<T-m) 

Similar in spirit calculation yields 

B. = ^^^ E su{a{t-t'))st.,+u{-mrt''\t)- (A.13) 



('m^ Scalar products. Now we are ready to express the scalar products in Eq. (A. 4) 
in terms of the r function Eq. (3.1). Having in mind the Heine formula Eq. (A. 7), we 
rewrite the scalar product in the l.h.s. of Eq. (A. 4) as 

P]"' (t; A)\ = 7^-^ / TT {d^j ^n{\o) e<*'^A ■ Aj{X) 

i 

Due to the identity 

I 

AKA)n(Am-A,) = A^+i(A), (A.15) 

the scalar product is further reduced to 

^^i— / TT fdA, r„(A,)e''(*;^^)) • A2 (A) , . (A. 16) 

l\rf-\i)k^^^]t^ ' } ^ n-=i(Am-A,) 

Symmetrising the integrand, 

n ■=l(A.+ l - A,) ^ ^ S n -tl, ,-^a(Aa - A,) ^"^'^^^ 
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and employing the remarkable relation (taken at n = ^ + 1) § 

Aa-A, I ls,(p„(A)), a>0 



a=l y j=l,j-5^a 

we deduce: 

f 0, a < 

For a > i (which is the only nontrivial case), this scalar product reduces to 

1 



A" 



By the same token. 



(Hi) The bilinear identity. Now we are in position to derive the bilinear identity. 
To this end wc substitute Eqs. (A.12), (A.13), (A.20) and (A.21) into Eq. (A.4). Up 

to the prefactor 

1 

its l.h.s. reads: 

oo 

j:A.{x^\pt\t;X)l 

<T=i 

CXD (T 

{{a - l)(t - t')) Sm-.+k {-[df]) s.-e {[dt]) rltt'""^(*)^i"""H*')- 

(T=^ A:=max(0, (T—m) 

Owing to the identity \\ 

■^fc(-[9t'])Ti™-")(t')=0 for k>m, (A.22) 
the latter reduces to 

oo oo 

^ ^ afe ((a - l)(t - t')) sm-a+k (-[a*']) s,.e {[dt]) ri^t'""^(*)^i"^-"^(*')- 

(T— £ fc—inax(0, CT — m) 

Interchange of the two series brings 

oo k-\-m—£ 

Sk {{a - l){t - t')) Yl Sk+m-e-a {-[dt']) {[dt]) T'l%'--\t)T^X--){t'). 

A;— max(0,£— m) ct— 

§ Equation (A.18) is essentially Eq. (3.16) written in terms of the Schur polynomials. 
II Equation (A.22) follows from Eq. (A. 11) and the expansion 

oo 

= Y^X^Sk{-PmW). 
k=0 

Indeed, since the l.h.s. is the polynomial in x of the degree m, the Schur polynomials (— p^(A)) in 
the r.h.s. must nullify for k > m. 



Y\ (1 — a;Aj) = exp 



3 = 1 



J = l 
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Making use of yet another identity J 

k 

Sk{t -t') = Y^ Sk-a{-t') s„{t), (A.23) 



CT=0 

we derive: 



oo 



t'=t 



5^ Sfe {{a - l){t - t')) su+m-e {[dt] - [df]) rS'""^(*)^i"-"n*')- (A.24) 

/i:— max(0,£— m) 

Setting the vectors t and t' to 

(t,t') 1-^ {t + x,t- x), 

(see Eq. (3.28)), and applying the Property 3 of Appendix B [Eq. (B.7)], we rewrite 
Eq. (A.24) in terms of Hirota operators: 

oo 

e^"-''^ E - ^I'+m-e {[D]) ri|t'""^(t) o rt-'-Ht). (A.25) 

fc— max(0,^— m) 

Finally, replacing n with n = m — s, we reproduce the l.h.s. of Eq. (3.34). The r.h.s. of 

Eq. (3.34) can be reproduced along the same lines starting with the r.h.s. of Eq. (A. 4). 
This ends the alternative proof of the bilinear identity in Hirota form. 

B. Hirota Operators 

Definition. Let f{t) = f{ti,t2,---) and g{t) = g{ti,t2,---) be difFcrcntiable 
multivariate functions. The Hirota derivative is a bilinear differential operator, 
Dkf{t)og{t), defined by 

D,f{t)ogit) ^^f{t + x)g{t-x) =(±-±.) f{t)g{t') 
oxk x=o \otk at).) 

Property 1. Let ?(£)) be a multivariate function defined on an infinite set of 
differential operators JD = {Di,D2, ■ ■■) such that 

?(£>) m o g{t) = T(a^) fit + ^)g{t - 4) 1^^^ (B.2) 

with = {d/d^i,d/d^2, ■ ■ •)• Then it holds: 

?(£>) m o g{t) = Vi-D) g{t) o f{t). (B.3) 

Property 2a. Setting g{t) = f{t) in Eq. (B.3), one observes: 

n^r/W°/(t) = for Y,ak=odd. (B.4) 

k k 

Property 2b. It follows from Eq. (B.l) that 

DkDef{t)o fit) = 2f\t) log/(t) (B.5) 

X It readily follows from the equality 

Taylor-expanded around x = and the definition Eq. (2.11) of the Schur polynomials. 
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and 

D.DeDMit)om = 2m g^^g^^^ log/(t) 



+4/'(*) 



92 \ / 



.2 \ / 92 



+ (sS;: (sk: '"^•f"' 



(B.6) 



Property 3. Let (x • D) denote the scalar product (x • D) = J^'h^i^kDk- The 
exponential identity holds: 

e(-^) /(t)o 5(t) = /(t + x) g{t - x). (B.7) 

For a more exhaustive list of the properties of Hirota differential operators, the reader 
is referred to Hirota's book (Hirota 2004). 

C. Jacobi Unitary Ensemble (JUE) 

The correlation function of characteristic polynomials in JUE is defined by the formula 

nn|p(^;«) = ^ / ilidXi{i-\^Y{i + \irf{{^.-XiYA-i^l{x) (c.i) 

" ^(-1,1)" J=l V Q,= l / 



where 

Tij + i)r(j + /x)r(j + 

n Tij + n + fi + u) 

is the normalisation constant. It is assumed that both /x > — 1 and 1/ > — 1. The 
associated r function equals 

tJ(^,k;*) = 1 /' n [dX, (1 - X.ni + Xjye^'-'-'^^ [] (,„ - A,)- ) ■ Al{X). 

(C.3) 

The superscript J standing for JUE will be omitted from now on. 

Although the very same technology is at work for a nonperturbative calculation 
of the JUE correlation function Eq. (C.I), its treatment becomes significantly more 
cumbersome. For this reason, only the final results of the calculations will be presented 
below. 

C.I. Virasoro constraints 

In the notation of Section 3, the definition Eq. (C.I) implies that || 

/(A) = 1-A2 ^ afe = 4,o-4,2, (C.4) 

g{X) = {fi - i^) + {H + v) X H> 6fc = (/i - i^)4,o + (At + z^)4,i, (C.5) 
D = (-l,+l) ^ dim(c') = 0. (C.6) 

II Notice that dim(c') = follows from Eq. (3.70) in which Zq = {-1, +1}. 
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This brings the Virasoro constraints Eqs. (2.15) (2.21) for the associated r function 
Eq. (C.3) in the form 

1 g 9+2 g 



m=0 



d 



«-™ m=0 

d 



dt, 



dt 



g+2 



9tq+2-m 

Tn{';,n,;t) = 'bqTn{';,n;t), (C.7) 



where the operator i)q is defined through operators Eq. (4.8) as 



a=l 



_9_ 
9?a 



and 'Cg(t) is the Virasoro operator given by Eq. (2.21). 
The three lowest Virasoro constraints for g = — 1, 5 = 0, and g = +1 read: 
[ 12^^^ ~ IZ^^^mTT: - i2n + ij, + u + k)—] logr„(^, k; t) 



0=2 



dti 



(C.8) 



+nti - n\pi-v+ 'di{(;,K)] = D_i logr„(<j, k; t), (C.9) 



EOO oo ^ 



-{2n + iJb + v + k)—\ logr„(<j,K;;t) 



logr„(<?,K.;t) +n[n + K- t?2(<i,'«)] = i>o log r„(<?,K.;t), (C.IO) 



and 



E ^'^^ 7)^ - E + [2n + - i?2 (^, /^)] ^ - 2- 



dti dtidt2 



d d \ 

-\pL-v+ ■!?i(<j,k;)] ^ - (2n + /U + z/ + K) — I logT„(<j,K;t) 

+n[i9i(^,K) - i93(^,'«)] 



5 



dti 



d 



2 ( ^log'^"(^''*;*) ) ( ^logr„(^,/c;t) ) = Di logr„(<i, k; t). 



(C.ll) 



C.2. Toda Lattice equation 

Projecting the first Toda Lattice equation Eq. (4.15) onto the hyperplane t = with the 

help of the first [Eq. (C.9)] and second [Eq. (C.IO)] Virasoro constraints, one derives: 

— J 
TLi : 

(2n + M + 1/ + K)^(i)?.i +i)o)-{tt + i^ + K)[n-u+ i?i(^:, k)] f)-i logn„|p(^; k) 
+ fogn„|p(^;K))' 

+n{n + ^ + u + k) {2n + h + v + k)"^ - [fi - v + 'di{<;, k)^^ 
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{2n + n + iy + nf [(2n + ji + v + nf - l] 
(2n + M + vf [(2n + /i + vf - 1] 

n„+i|p(^; /«) n„_i|p(^; k) 



X 4n(n + /u)(n + z^)(n + + z^) ■ 



(C.12) 



Notice the structural similarity between TL^ and the first Toda Lattice equation for 
the CyUE [Eq. (D.12)]. 

C.3. KP equation and Painleve VI 

Projecting Eq. (4.57) onto t = with the help of all three Virasoro constraints 
Eqs. (C.9) - (C.ll), one derives: 



KP,: 



logn„|p(^;K) 



+ (l)^ilogn„|j,(^;K 
+ (l)_ilogn„|p(^;K)) 



+ 4230 



logn„|p(^;K) 



4DoD_i +2Di - 2D_i 



logn„|p(^;K)=7. (C.13) 



Here, Q;[j,fe], /3j and 7 are the short-hand notation for the following functions: 



"[-1,-1] 

"[0,-1] 
"[0,0] 

"[1,-1] 

^0 

/32 



= 4n2 + 2 - [//, - ;/ + ?9i(<;, k)]^ -4{h + u + k)[k- '!?2(^, k) - n] , 
= - 2(/i + z/ + k) - + i?i(<j, k)] , 

= 3 (2n + + + k)^ , 



= 2 



1 - 2 (2n + /X + 1/ + k) 



= (/Lt + + k) [/i - J/ + Z?3(^, k)] + [/i - + k)] [k - t?2(^,K)] , 

= - 2 [(2n + At + J/ + + (/i + 1/ + k) [k - t?2(^, k)]] , 
= - + t^ + k), 
= 2 (2n + + 1/ + k)^ , 

= 2n(n + /It + 1/ + k) [k - '!?2(<?, k)]^ + (/x + + k) [k - i?2(<i, «)] 

+ [/X - + k)] k) - z?3(^, n)] 



Remark. For p= 1, the equation KP^ simplifies. Indeed, introducing the function 



d 



<P(0 = (^^-l)^logn„(^;«;)+^ ^ bibj-bij^bj' 



(C.14) 



2<j<j<4 



j=2 



where 6j's are given by 



61 = +^{fJ'-i^), 



62 



2(m + «^), 



^3= +-(/x + zv + 2n). 



&4 



-(M + z/ + 2n + 2«;), 



(C.15) 
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one observes that Eq. (C.13) transforms to 

(^2 _ 1)2^/,/ ^ 2c(^2 - 1) if" + 6(^2 - 1) {ip'f - S><;W + 2<^' 

- Aui Lp' - 4I/4C - 2z/2 = (C.16) 

with the parameters 

4 

^1=13^'. ^2=^6^62, Z.3= ^ 62^2^2^ J/4=6iM3&4. (C.17) 
j=l j<j i<j<k 

Equation (C.16) can equivalently be written as 

- 4z/i {ip'f - Su4 -ip]- Au2ip' - Avs = 0. (C.18) 
The boundary condition at infinity reads: 

^--U + v + 2n)\{l + 0{<;-^)). (C.19) 

It is easy to convince yourself that Eq. (C.18) can be reduced to the sixth Painleve 
transcendent. Introducing the new function 

h{t) = \^{,) , (C.20) 

Z <;=2t — 1 

one straightforwardly arrives at the a form of the sixth Painleve transcendent (Forrester 
and Witte 2004): 

4 

h' [tit - 1) h"f + [h' - {2t - i)h') + hb^hbif = n('^' + ^')' c^-^i) 

see Appendix E for more details. 

D. Cauchy Unitary Ensemble (CyUE) 

For the Cauchy Unitary Ensemble (CyUE), the correlation function of characteristic 
polynomials is defined by the formula 

where 

cy ^ (27r)" fr r(j + l)r(2z.+ l-n-i) 

2"(2^-")jy r2(z/ + l-j) ^ 

is the normalisation constant. It is assumed that the parameter v is real and 
u > n + {n — l)/2, where n — tip k. The associated r function equals 

r„«^(<,,K;t) = 1 y^^ n (tttaif n(^« - ^^^^^ j • ^-w- 

The superscript Cy standing for CyUE will further be omitted. Similarly to the JUE 
case, involved character of the calculations makes us present only the main results. 
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D.l. Virasoro constraints 

In the notation of Section 3, the definition Eq. (D.I) imphes that 

/(A) = I + A2 ^ ak = Sk,o + Sk,2, (D.4) 

g{\) = 2uX ^ bk = 2uSk,i, (D.5) 

D = M ^ dim(c') = 0. (D.6) 

This brings the Virasoro constraints Eqs. (2.15) - (2.21) for the associated r function 
Eq. (D.3) in the form 

d _ ^ d 



d 



-2u- 



Tn(s, K; t) = Tlq Tn{<;, K; t), (D.7) 



dtg + 2 

where is defined through operators Eq. (4.8) as 

■D, = 'B, + ■B,+2 = E(l + ^r^^ (D.8) 

a—1 

(not to be confused with the operator Tig defined by Eq. (C.8) of the previous subsec- 
tion) and £g(t) is the Virasoro operator given by Eq. (2.21). 

The three lowest Virasoro constraints for q = —I, q = 0, and q = +1 read: 
E + E 5^;;^ - - 2n - /.) ^ log r„ n; t) 

j=2 ■> j=l ■'^ / 

+nti + ndi{(;,K) = logr„(^, k; t), (D.9) 



,J=1 J=l -1^ >- 



d \ ( d 

-{2u-2n-K)— I logr„(^,K;;t) + ( — logr„(<j, k; t) 

W+n[K+ ■&2{';,i^)] = i>o logT„(<?,K;t), (D.IO) 

(°° Q °° Q g 

+ ^,(., K) A + (. + 2n - 2.) A + 2^) logr„(., t) 
+n[t?i(^,K)+ t?3(^,'«)] 

+2(^^logT„(^,K;t)^ (^J-logT„(<i,K;t)^ = logT„(<i, k; t). (D.II) 
£1.^. Tbrfa Lattice equation 

Projecting the first Toda Lattice equation Eq. (4.15) onto the hyperplane t = with the 
help of the first [Eq. (D.9)] and second [Eq. (D.IO)] Virasoro constraints, one derives: 

Cy 

TL/: 
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{2u-2n-KyiVi^-^o)-i2iy-K)M<i,i^)^-i\ log n„|p(<r; k) 

= (2zy - 2n - Kf[{2i/ ~2n- nf - 1] 

n{2v-n) n„_i|p(^;/'i)n„+i|p(^;K) 

^ 4(n - - 1 ny^T^^ • 



(D.12) 



Notice the structural similarity between TL^ and the first Toda Lattice equation for 
the JUE [Eq. (C.12)]. 

D.3. KP equation and Painleve VI 

Projecting Eq. (4.57) onto t = with the help of all three Virasoro constraints 
Eqs. (D.9) (D.ll), one derives: 



KP 



Cy 



+ (i)^ilogn„|p(^;«)) 

- (^I)_ilogn„|j,(^;K) 



logn„|p(<?;K;) 

logn„|p(^;/«) = 7. (D.13) 



4DoI'-i - 2Di - 2'D_i 
Here, ct[j,fc], and 7 are the short-hand notation for the following functions: 
«[_!_!] = - 4n^ - 2 - - 4(2:/ -k)[k+ ^2{<;, K)-n], 

a[o-i] = 2(2i/-k;) 

3 (2z/ - 2n - k)^ , 



"[0,0] 

^0 
^1 

7 



= 2 



1 - 2 (2zy - 2n - K,y 

= {2u - k) i?3(<?, «) + «) [k + i?2(s:, k)] , 
= - 2 [(2z/ - 2n - - {2v -«;)[«+ i?2(<?, «)]] , 

= - di{c;,K){2v- k), 
= -2(2z/-2n-K)^, 

= 2n{2v -n-K) - [k + ^?2(<i, n)f + (2;/ - k) [k + a?2('?, «)] 

+ l?i(^,K)[t?i(^,K)+ l?3('?,«)] 

Cy 

Remark. For p = 1, the equation KPj^ can be simplified. Indeed, introducing the 
function 

d ^ 
(^(c) = (l+^')^logn„(^;K)+^ ^ 5i6,-6i^6„ (D.14) 



2<j<j<4 



j=2 



where 6j's are given by 



61 =0, 

h2 = v-n, 

63 = - I/, 

hi = K — u -\- n, 



(D.15) 
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one observes that Eq. (D.13) takes the form of the Chazy I equation 
(1 + + 2c(l + c') + 6(1 + {ip'f - 8c V'V' + 2v^2 

+4z/i + 2!/2 = (D.16) 

with the parameters 

0=1 i<j i<j<k 

Equation (D.16) can cquivalently be written as 

[(1 + ?2)<^"] 2 + + V(?(^' - + Aut + Au2<f' + 4!/3 = 0. (D.18) 

The boundary condition at infinity reads: 

<p(0 ^-{n-v)\{l + 0{^-^)). (D.19) 
Introducing the new function 

one straightforwardly verifies that it satisfies the equation 

4 

[t{t - 1) h"f + h' [2h - {2t - l)h'f = \{{h' + h]). (D.21) 

This coincides with the a form of the sixth Painleve transcendent [Eq. (C.21)] with 
6i =0. 

E. Chazy I Equation as a Master Painleve Equation 

This Appendix, based on Chapter 6 of Forrester (2010) and Section 4.3 of Adler and 
van Moerbeke (2001), collects very basic facts on six Painleve transcendents and a 
closely related differential equation belonging to Chazy I class (Cosgrove and Scoufis 
1993). 

E.l. Painleve property and six Painleve transcendents 

Painleve transcendents are second order differential equations of the form 

V" = F{t,y,v'), (E.l) 

for which all movable singularities of y{t) are limited to poles, given is a rational 
function in all its arguments. Note that this requirement, known as the Painleve 
property^ does not rule out the existence of immovable (i.e. fixed) essential singularities. 
We remind that a singularity is called movable if its location depends on one or more 
integration constants. 

Painleve (1900, 1902) and Gambler (1910) have shown that the requirement that all 
movable singularities arc restricted to poles leads to 50 types of equations, six of which 
cannot further be reduced to either (i) linear second order differential equations or (ii) 
the differential equation for the Weierstrass ^-function, 

{y'f = 42/3 - g^y - g^, (E.2) 

{92,93 are constants) or (iii) the Riccati equation 

y' ^a{t)y^ + b{t)y + c{t), (E.3) 

where a{t), b{t) and c{t) are analytic functions oit. (The two later equations represent 
an irreducible class of first order differential equations of the form P{t,y,y') =0 with 
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P being a polynomial in y' , y with coefficients meromorphic in t, such that y{t) is free 
from movable essential singularities.) 



Labeled Pi to Pyi, these six equations are explicitly given by 



Pi ^ 
Pii 

Piii ^ 

PiV: ^ 



VI 



y" = 6y2 + 1, 

y" = 2y^ + ty + a, 

y" = -{y'f-jy' + ay' + Upy' + ^) + -, 

y t t ^ ' y 

y" = Yy ^y'f + ^ + 4i + 2(i^ - «)y + ^, 



Pv ^ y" = 



1 



2y \-y 



{yf 



(y-i)' 



ay ■ 



y 



ry Sy{y + 1) 
^ t ^ y-1 ■ 



-(- 

2 \y 



y - 1 y-t 



{y'V 



1 



1 



t-l y~t 



, y{y-l){y-t) ( Bt 5t{t - 1) 

H —. — — — I a + — + 7 7^ + 



(E.4) 
(E.5) 

(E.6) 
(E.7) 



(E.8) 



(E.9) 



t2{t-iY V y^ [y-if [y-ty 

Here, a, p, 7 and 5 denote complex constants. The equations Pj, Pu, and Pjy have 
essential singularities at the point 00, the equations Pu and Py have critical points 
and 00, whilst the Pvi equation has critical points at 0, 1, and 00. 



E.2. Hamiltonian formulation of of Painleve transcendents and their 
Jimbo-Miwa-Okamoto a forms 

Each Pj of Pi - Pyi is equivalent to a Hamiltonian system {Q,P,Hj} 

dHj 



Q = + 



dP '■ 
dQ ' 



(E.IO) 



where Q and P denote derivatives with respect to t, and the Hamiltonian _ffj = 
Hj{P, Q, t} is a polynomial or rational function. The equation in canonical coordinate 
Q, obtained from Eq. (E.IO) after eliminating canonical momentum P, is appropriate 
Pj equation. Existence of the Hamiltonian formulation of Painleve equations can be 
traced back to absence of movable branch points in Pj (Malmquist 1922). 

Explicit forms of the Hamiltonians Hj, as well as the differential equations satisfied 
by them, can be found in original papers of Okamoto (1980a, 1980b), see also Noumi 
(2004) and Chapter 6 of Forrester (2010). In the random-matrix-theory literature, the 
differential equations for H] more often appear in the so-called Jimbo-Miwa-Okamoto 
(T-form: 



(an' + 4(a0"-2(ai-to;)-0. 



(7 Pi 

crPii (fTn) 



4a(T 



{a'uf - ta[i + Gil 



- 6^ = 0, 



o-Piii [tcTiiif + cTiii(4a-{ii - l)((7iii - toiii) 



2 2 



(E.ll) 
(E.12) 

(E.13) 
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aPiv ^ (a{'v)' - 4(aiv - ta'.yf + Aa'^y {a',y + bj) = 0, (E.14) 

4 2 4 

aPv ^ {ta'^f - [ay - ta'y + 2{a'yf + ^ 6,] + 4 (^V + = 0, (E.15) 

i=i i=i 

4 2 



i=i 



(E.16) 



-nKi+^')=o. 

These follow from equations for the polynomial Hamiltonians after appropriate linear 
transformations of Hj (Forrester 2010). 

E.3. Chazy classes and a master Painleve equation 

Differential equation belonging to Chazy classes (Chazy 1911) is the third order 
differential equation of the form 

y'" = F{t,y,y',y") (E.17) 

with F being a rational function in y, y' , and y" and locally analytic in t, whose general 
solution is free of movable branch points. In his classification, Chazy discovered 13 
irreducible classes (reviewed in Section 6 of Cosgrove (2000)), of which the class Chazy 
I is given by 

„, P', 6,0 4P' , P" 2 4Q , 2Q' 2P „ 

with arbitrary polynomials P{t), Q{t) and R{t) of third, second and first degree, 
respectively. 

Cosgrove and Scoufis (1993) have proven that Chazy I equation (E.18) admits the 
first integral, which is of second order in y and quadratic in y". 



iy"Y + p, 



{Py'^ + Qy' + R)y'- {P'y'^ + Q'y' + R') y 



■i(PV + Q")2/'-^^V + c 



0, 



(E.19) 



c is the integration constant. They also show that this equation is the 'master Painleve 
equation' because it unifies all of the six Painleve transcendents into a single equation. 

Making use of Eqs. (E.18) and (E.19) in conjunction with Eqs. (E.ll) (E.16), one 
readily verifies that uPj can be brought to the Chazy I form (CPj): 

CPi ^ < + 6(aJ)2 + t = 0, 

CPii ^ < + 6(aii)2 - 4^C7^I + 2cTii = 0, 

2 ^ 

CPiu ^ tVf/i + tai'ii - 6t {a'uif + Aama^u + + H "^"i ~ 2 "^"^ " °' ^^'^^^ 



(E.20) 
(E.21) 



CPiv 
CPv 



i=i 



a[y + UaiY + 2 JJ 6j = 0, 



(E.23) 
(E.24) 
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+ 



3=1 l<i<j<4 J j=l 

CPvi ^ [t{t - 1)] V^'i + t{t - l)(2t - 1) a^^i + 6t{t - 1) {a'yif 



(E.25) 



4(2t - 1) aviaVi + a^i ^ 6^ + 2a^, 



+ 



* E ^^''1+ E ^^^^^ 

l<i<j<4 l<i<j<fc<4 



= 0. 



(E.26) 



F. Functions of Parabolic Cylinder and Two Related Integrals 

In this Appendix, we collect some useful formulae related to the functions of parabolic 
cylinder and also treat two related integrals Eqs. (5.27) and (5.75) encountered in Sec- 
tion 5. 



(i) Integral representations, differential equation, and a Wronskian. — The functions 
of parabolic cylinder Dp{w) and D-p-i{iw) admit integral representations 

1 



VTT 



-x^ /2—ixw 



and are general, linear independent solutions to the equation 



9\tp > -1, 
Olep > -1, 



dw'^ 



The corresponding Wronskian is given by 
V^^[D_p_i{w),Dp{iw)] = 



(F.l) 
(F.2) 

(F.3) 
(F.4) 



(ii) Functional and recurrence identitites. — The following functional identity holds: 

(F.5) 



Dp{iw) = ^^^^J^ [(-i)f£»_p_i(w;) (-«;)] . 



277 

In addition to the three term recurrence relation, 

Dp+i{w) - wDp{w) +pDp_i{w) = 0, 

there exist two differential recurrence relations: 

w 

D'p{w) + -Dp{w)-pDp_i{w) = 0, 
D'piw) - -Dp{w) + Dp+i{w) = 0. 



(F.6) 

(F.7) 
(F.8) 



(Hi) Asymptotic expansions. — Throughout the paper we make use of the large-iz; asymp- 
totic expansions: 



• For |argw| < 37r/4: 



1 + 



(F.9) 
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For 7r/4 < argw < 57r/4: 



Dp{w) ~ w^e 



,Pp-™ /4 



l + 0( — 



n-p) 



1 + 01 — 



(F.IO) 



For — 57r/4 < argw < — 7r/4: 



1 + 



r(-p) 



1 + 



(F.ll) 



(iv) Integrals Eq. (5.27) and Eq. (5.75). -To calculate the integrals Eq. (5.27) and 
Eq. (5.75) encountered in Section 5, we first formulate the Lemma. 

Lemma. Let Ui{t) and U2{t) be linearly independent functions whose Wronskian 
is a constant: 

^([1*1,^2] = U1U2 — u[u2 = Wo ^0- (F-12) 

Then, 

dt a-iui{t) + a2U2{t) 



I 



(F.13) 



(ciUl(t) + C2U2(t))^ ClUl(t) + C2U2{t) ' 

where ai and a2 are the constants satisfying the relation 

cia2 — C20L\ = —. (F.14) 

Wo 

Proof. Differentiate both sides of Eq. (F.13) and make use of Eq. (F.14). 

Remark 1. The integral Eq. (5.27) follows from the Lemma upon the choice 

ui{t) = D^N-i{t), U2{t)=DN{it). (F.15) 
Since wo , we conclude that 

dt aiD-fq-i{t) + a2Dff{it) 



{cxD_M-l{t) + C2DN{it)f CiD_N-l{t)+C2DN{it)' ^^'^^^ 

where 

ciQ!2 - C2ai = {-i)^. (F.17) 



Remark 2. The integral Eq. (5.75) follows from the Lemma upon the choice 

ui{t) = DN-i(f), U2{t)=D_N{it). (F.18) 
Since wq = {—i)^, we conclude that 
f dt aiDN-ijt) + a2D-N{it) . 

J {ciDN-l{t)+C2D_N{it)f ~ CiDN-l{t) + C2D_N{it)' 

where 

cia2 - C201X = i^. (F.20) 
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